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NATIONAL DEVELOPMENTS 


IMPLEMENTATION OF TECHN(CAL CONTRACTING SYSTEM DISCUSSED 


Beijing KEYAN GUANLI [SCIENCE RESEARCH MANAGEMENT] in Chinese No 3, 1986 
pp 28-32 


[Article by Chen Changhui [7115 2490 1741] of the Science and Technology 
Bureau, Ministry of Aerospace (Industry: "On Problems Associated with the 
Implementation of Economic Responsibility System for Technical Programs at 
Research Institutes") 


[Text] The reform of the funding system to research institutes is an 
important measure in the reform of the S&T system. In order to adapt to this 
refora, research institutes aust implement various technical economic 
responsibility system centered around contracting bdased on  acientific 
@anagement methods. The following is a discussion of several problems 
associated with the iaplementation of this technical economic responsibility 
Systex by research institutes. 


I. Necessity of [Implementing Technical Economic Responsibility System 


In the past, the Chinese government paid for all research expenses. There are 
serious drawbacks associated with the distribution of research money since 
everyone is getting funded from the same pool of funds. 


1. There is no economic pressure at research institutes. The assignments are 
handed down from above and the funding also comes from above. The 
accomplishments are also reported to superior organizations. Thus, technology 
is detached from economics and scientific research is separated from 
production. This affects the wide application of research accomplishments. 


2. There is a lack of scientific basis and rigorous proof in the distribution 
of funds. Most of the time, the goney is split based on historic levels. 
Thus, the research assignments are not in sync with fund distribution, 
Projects with higher allocations often spend the extra money on other items. 
Projects with insufficient funding cannot keep with the schedule. 


3. Research funding is managed out of phase with spending. The technical 
personnel only spend and do not manage. Whatever is spent is getting 
reimbursed. Research funding is not rigorously controlled. However, 














research expenses and materials consumed are somewhat flexible in nature. 
Waste will occur if the technical staff does not consider the economics. 


In order to overcome these problems, the funding system must be changed. 
Externally, we must implement a technical contract system. Internally, we 
must have a techrical economic responsibility system. 


This reform was implemented in several experimental sites for a few years. 
The experience shows that the management of research funding is improved and 
some research money is saved. The research people have a stronger sense of 
responsibility and the efficiency of research is improved. For instance, 
after implementing the responsibility system for the development of a certain 
product in our ministry, the research staff followed a strict budget to save 
expenses. In addition, the organization was managed based on system 
engineering principles. The result is that the development time is reduced by 
two thirds as compared to that of a similar project prior to the 
implementation of the responsibility system. The cost is saved by one half as 
well. The implementation of this system also encourages the research 
institutes to be oriented toward the economy and the society. It will break 
them away from department ownership and locality ownership, end will alter the 
various dependencies under the present structure. It will promote the Linkage 
between research and production to allow technical accomplishments to be 
rapidly applied to manufacturing so that they can be converted to profits and 
social benefits. In the mean time, new sources for research funding are 
created. For example, a certain research institute in our ministry had a pure 
annual profit of approximately 2,000,000 yuan in 1983 and 1984. Due to the 
implementation of the reform, the 1985 profit is expected to be over 5,000,000 
yuan. The institute created a research and development fund and a collective 
welfare fund from the profit. The government has reduced its support by 15 
percent. 


II. Methods for Implementing a Technical Economic Responsibility System 


To implement a technical economic responsibility system and external contract 
System, we must bring economic and management concepts into the domain of 
scientific research to directly link research assignments with expenses. The 
technical and economic responsibilities of the research unit must be clearly 
identified so that responsibility, power and benefit are combined. 


The ministry issues the research assignments to the institutes. The ministry 
Signs a primary contract with the lead organization which then signs various 
lateral subcontracts with other collaborating units (or subcontractors of 
various systems). The contract aust specify the tasks, accomplishments, 
schedules, costs and bonuses. Under the premise that the completion of ali 
contracts within the ministry must be ensured, the institute can sign lateral 
contracts with outside organizations to face the economy and the society. 
Within the institute, all contracting responsibilities are assigned to various 
levels and individuals. The following is a discussion on the methods employed 
to implement the economic responsibility system in the institute. 
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1. General Contract Method 


Different contracting methods should be applied to different research 
laboratories and shops. The situation is introduced in the following: 


a. Major laboratories and machine shops, because of their ability to undertake 
research and production assignments of specific directives and to accept 
external contracted work, have their own incomes. They should be treated as 
independent accounting units, responsible for their own profits. The 
contracts should primarily contain targets specifying the research and 
production items as well as incomes. The research funding is determined by 
subtracting a specific proportion of the contracted amount for a risk fund and 
for management fee. At the end of the year, or upon completion of the 
project, a portion is set aside as bonuses based on the status of the 
assignment and on tne net income. There are two ways to set the bonuses. One 
is to assign different percentages to net incomes derived from contracted 
targets and additional objectives, respectively. This is done either to 
encourage bringing in more work, to receive more income, or to ensure the 
completion of research projects in the plan. Rach unit can decide its own 
priority based on its own situation. The other way is to bundle planned anc 
unplanned projects together to arrive at a goal. Only one percentage is set 
for bonus payoff. This is a more convenient way to manage projects. in 
addition, it can prevent laboratories and shops from diverting costs ineurred 
in contracted work of their own to projects planned by the institute. 


b. A total duty responsibility system can be implemented with service-oriented 
laboratories and shops such as informition research offices, instrument 
laboratories, chemistry laboratories and machine shops. (Their assignments are 
Mostly service-oriented in support of major laboratories and shops in the 
institute. They can also accept outside work to receive technical service and 
service income.) The personnel, positions, and resposibilities are fixed. 
They receive compensation for some of the services rendered. On the basis of 
overall evaluation, bonus is given based on the completion status of the 
contracted tasks. A portion of the income from external services rendered is 
set aside for the employees as well. 


2. Contract Method Based on the Implementation of Economic Objective 
Management 


Because of the need to gradually reduce government support to achieve economic 
independence, research institutes in the technology development agode can 
implement the economic objective management system. The net income of the 
institute should be determined based on the amount of operating expenses to be 
saved each year. Based on factors such as the total wages of the people 
working on the first Line of research and production, fixed capital, oecupted 
area and liquid assets, the economic goal of the entire institute is broken 
down level by level to determine the amount to be undertaken by each unit, 
including directive research and production assignments. This mothol can 
stimulate the research laboratories and machine shops to actively tap into 
potential resources, 3ave expenses, and add now projects to make the institute 
economically independent. Nevertheless, too much stress on economle 
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objectives can make the institute go astray from 
research. 


3- Contracting Single Project 


Some important researc.: projects which involve several departments can 
contracted to the task groups. By implementing the system of fixing the 
assignment, accomplishments, schedule, expenses and bonus in the contract, 
these projects can be directly managed by the institute. if other 
labdoratories or shops are involved in some tasxs, they can sign subcontracts 
with the task teas. The leaders of such teams must be chosen by the 
institute. When appropriate, we can also try a bidding method. 

IIil. Organization Management with Implementation of Technical Econom! 
Responsibility Systea 


The implementation of a technical economic responsibility system with laters 

and serial contracts for research projects will make {t easier to organize 
research based on systems engineering methods. It will also facilitate the 
coordination of various tasks. However, we should not neglect management 
because we implement the contract system. m the contrary, because economic 
ang Management concepts are brought into the domain of research, management 
will have to cover a wider range. The work becomes more complicated whic? 
requires better management. Specifically, we have to pay attention to the 
following points: 


First, in the past, our primary concern was focused on technical feasinility 
and applicability and was not concerned with economics. To implement the 
technical economic contract system is to combine both for proof of technology 
and economics in order to ensure the accuracy of the research budget and the 
economic benefit of the project. 


Second, there are more sources for projects after we make research face the 
society which adds a new dimension to organization management. The 
utilization and ba'‘ance of research resources becomes more complex. 


Third, as we implement the system to combine responsibility, authority and 
economic interest, %o link technical accomplishments to economic gains, the 
personal welfare of the technical staff is affected. When the work is 
contracted, it will require tighter management and more rigorous § inspection 
and review. 


Fourth, due to the exploratory and uncertain nature of research, no matter how 
detailed the contracting work is done, there will be many unforeseen problems 
which require immediate coordination by the management. 


Fifth, after implementing the contract system, each department is more 
concerned about its own interests and tends to neglect the overail picture. 
For instance, departments with general processing and testing equipment are 
happy to accept. external work for additional income. They are not 
enthusiastic about assignments within the institute. Some departments will 
overload the equipment for more money. Some task groups would rather 0 








outside because they feel the institute charges too guch. Thus, the institut« 
loses some funding. In order to save money, some task teams do not do enoug 
experiments. Therefore, we must perfect the structure of the institute wit! 
specific reference to the problems emerging after the implementation of tix 
sontract system so that internally it is balanced and coordinated. Specia 
favorable internal processing and testing rates must be established to promot: 
collaboration and unity among all departments so that the institute is stil. 
one entity. In addition, management of technical quality sust be 
strengthened. 


implementation of the tecnnical economic responsibility system involve: 
everyone. It requires the strengthening of three functional departments-- 
planning, finance and maintenance. Furthermore, we sust establish ar 
effective technical control and command systen. 


The research planning department must handle the relations between what’ 
needed and what's feasible, overall and local, focal points and genera 
issues, progress and quality, and long term and short term well. It plays 3 
leading role in assigning tasks, issuing contracts and rewardin 
accomplishments. The finance department sust regulate, gonitor and serve th 
institute to organize and manage the cost accounting work. in order to matc 
planning with economic management to allow the proper management an 
utilzation of research funding, we need dedicated organizations and personne 
to manage research funds. In general, we may set up a technical economic 
office under the planning department. The maintenance department mus' 
properly manage the ir 3truments and equipment to meet the requirements of th 
new system. 


In order to make sure that technical and quality-related management work su: 
as development planning, project demonstration, and technical coordination and 
review is done properly to ensure the quality of research, we cannot weake: 
the technical command functions headed by the chief engineer and chie: 
designer in implementing the contract systen. Moreover, it requires a clos 
coordination with the administrative management systen. 


The overall research and production dispatch system is an important part < 
the administrative command system. Its duty is to integrate th 
responsibility and authority of the department which ensures the resources 
needed in the research into the assignment to effectively utilize them to 
guarantee the success of the project. Amain dispatch office may bd» 
established under the research planning and management department. It may 
comprised of leaders from the planning department, maintenance departments, 
laboratories and shops to form a central command system centered arour 
research and production. The chief dispatcher, acting on behalf of the 
director of the institute, exercises the power to issue orders on regular or 
emergency basis to ensure the completion of the assignments. A temporary 
emergency order may be considered as a supplement to the original contract. 
Emergency orders should be approved by the director of the institute. 


In various stages of research, we must closely combine the technology wit! 
economics. In the proof of feasibility stage, we should fully demonstrat» 
both technical and economic feasibility by involving technical, planning anc 














economic talents to choose the optima) plan. in the design development stage, 
cost must be controlled by value engineering and cost-effect analysis methods. 
The allocation and utilization of research funding wust be broken down, 
tracked, accounted for and analyzed in order to reduce costs in various 
stages. In the accomplishment certification stage, not only technical review 
must be made, but also cconomic and social benefits must be evaluated. Only 
vy doing so, the technical economic responsibility system can be fully 
effective to improve the economic welfare of the institute. 


The implementation of the technica] economic responsibility system is an 
important step for a research institute which is changing from a pure research 
mode to a research and production management mode. Systematic reforms in the 
leadership system, organization structure, wage system, planning systen, 
internal relation structure and politica) ideology must be made to make the 
transformation effective. 


IV. Principles to Follow and Relations to Keep in the Reform 


When the institute implements reform, we must insist that the direction of 
research is the most important issue. We must «lways think based on the goals 
of achieving more accomplishmnts, training more professionals, and creating 
more social and economic benefits. We must place the assignments given by our 
government and our superiors first. In addition, we should pay attention to 
strengthening our political ideology and spiritual civilization. We must 
insist on the principle of distribution based on work. The following three 
relations must be properly handled: 


1. Proper Treatment of the Relation between Government, Collective Entity and 
Individual) in Implementation of Consolidation of Responsibility, Authority and 
intercst 


The core of the technical] economic responsibility system is “responsibility.” 
We must determine “authority” and “interest” based on “responsibility.” 


In the area of “responsibility,” in the reorganization within our ministry in 
recent years, a comprehensive technical post responsibility system has been 
set. up. It is a basis for implementing the technical economic responsibility 
system. However, it still must be perfected based on the new situation. In 
the area of “authority,” the Ministry issued a ruling in September 1984 to 
expand the autonomy of individual business units. More autonomy is given in 
sreas such as planning, personnel, finance and discipline. In terms of 
leadership system, most institutes employ a director responsibility system. 
In addition, tne authorities of research teams and shoos are also expanded. 
Some attempts were made to organize task teams by free choice. These measures 
also facilitated the implementation of the new contract systen. In the area 
of "profit" it is primarily done through keeping a fraction in the institute 
and distributing bonuses to individuals. Tne Chinese Government specifies the 
bonus exemption limit and the bonus taxation method. Each operating unit also 
has rules and regulations governing bonuses. In terms of bonus aliocation, 
the size of the bonus will certainly affect individual enthusiasm. 
Nevertheless, the key is to overcome the idea of even distribution. It should 
be allocated based on performance. The bonus gap should be widened 
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reasonat ly. Currently, bonus allocation is made based purely on profit. Due 
to the .mreasonable social value systez, different profit margins for each 
product, different environments, and other factors, this method may result in 
unequal rewards for equivalent amounts of work. Based on a survey, the 
difference is of the order of six- to eightfold. Such unreasonable gaps will 
affect the enthusiasm of some employees. in order to sinimize this gap and to 
@aintain relative internal balance, some institutes also collect a bonus tax 
from the basic departments based on the sethod implemented by the government 
to regulate the imbalances. 


When we handle the relation between our country, the institute and an 
individual, we must strengthen the political ideology work to rectify the 
Misconception that “goney is everything.” We gust combine the evaluation of 
the technical economic responsibility systems with the realization of spiritual 
civilization. We have to link the evaluation results with performance 
appraisal and promotion so that the consolidation of responsibility, authority 
and welfare is realized in a wider perspective. 


2. Careful Handling of the Relation between “Building Foundation” and 
"Catching up With the Norm” with Economic Benefit in Determining the Guiding 
Business Strategy 


After we opened up the institutes and expanded their authorities, there are 
two situations to be concerned about in the choice of tasks and in the guiding 
business strategy. One is that tasks related to “foundation building” and 
"achieving the nors," because of the degree of difficulty and length of 
Guration, are not economically attractive in the short term. WNobody wants to 
undertake such assignments. Everyone is interested in tasks with short cycl:, 
high profitability and of appropriate technical standards. In some cases, 
people abandoned research and got into business to gake soney. The other 
extreme is to be detached from current resources and blindly pursue leading 
high technology items by ignoring the economic benefits of the short tern, 
appropriate level, profitable projects. Both extremes gust be avoided. 


When the relation between “level” and “economic benefit" is treated, we must 
base on an analysis of the specific situation. Institutes with better 
foundation and resources, under the premise of facing economic construction 
and emphasizing economic benefit, sust strengthen their involvement in the 
exploration of new technology and focus on long range projects and technology 
reserves in order to gaintain their advantages. Institutes of weaker 
foundation and poorer resources gust be practical. They should realize that 
although the majority of those short term, appropriate level projects does not 
belong to high technology, yet they are small, highly specialized, highly 
profitable technologies needed to vitalize the economy. 


During the early stage, regardless of the foundation of the institute, the 
desire to improve the economic benefits is gore urgent. It is natural to do 
tne easy and profitable tasks first. This is not only a current need but also 
an unavoidable process in the transformation from a closed to an open type of 
research institute. During this process, we can adjust to the new management 
System, explore new avenues, get into new habits, and accumulate required 
experience. However, as the technical business grows, low level projects wiil 








be Jess competitive. Research institutes will be forced to raise the 
standards and improve technology reserves. 


3. Appropriate Treatment of HKelation between Fundamenta) Mznagement and 
Reform 


Sowe fundamental management is required to implement the technical economic 
responsibility systes. Through the reform, fundamental management can be 
further perfected. These two things are complementary. In reality, once we 
talk about “improve economic benefit" and "be flexible," we are anxious to 
make woncy. The existing management system and operating procedure which are 
still valid and current are being set aside. For instance, a project may be 
initiated without going through rigorous technical and economic proof of 
feasibility. In some cases, the quality control system is inadequate, leading 
to repeated work. In some cases, speed is the only concern and the 
accumulation of original technical and economic data is ignored. If these 
things are not rectified in time, the reorganization results over the past 
several years will be irrevocably lost. Scientific and economic laws must be 
obeyed for the reform to strengthen fundamental management work. 


There is much work in fundamental management. However, based on implementing 
the technical economic responsibility system, first we must perfect the 
technical post responsibility system, as well as other rules and regulations, 
according to the neods of reforms. In particular, we must incorporate contents 
in economic management and development to be implemented thoroughly. Second, 
we must strengthen the statistical work. We have to hire people, set up a 
system, perfect the research statistic indicator system, and establish the raw 
data forms and tables to create a complete statistical network. Then we must 
accurately calculate the cost. Through projection, planning, control, 
accounting, analysis snd review, the cost and economic benefit of a research 
product are reflected. We should promote to combine technical and economic 
consideraticns to save expenses and reduce cost. To do a good cost 
calculation, we need to do a thorough inventory, formulate an internal pricing 
structure and depreciation schedule, define contents covered by the audit, 
clarify the roles of each functional department in cost management and 
accounting, establish an audit system and internal billing system, and 
prepare to keep good original data. The Ministry has already formulated a 
cost managescnt and estimation method for research organi zations. It will be 
implemented on a trial basis on 1 January 1986. The cost estimation work at 
research units will be initiated by the implementation of this method to 
further perfect the technical economic rosponsibility systen. 
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APPLIED SCIENCES 


MICROCOMPUTER USE AS FURNACE CONTROLLER DESCRIBED 


Beijing DIANZI JISHU YINGYONG [APPLICATION OF ELECTRONIC TECHNOLOGY] in Chinese 
No 2, 25 Feb 86 pp 6-9 


[Article by Wang Dexi [3769 1795 3886] and Song Jian [1345 0256], Research 
Institute No 13, Ministry of the Electronics Industry: "The Application of 
Microcomputers in a High Temperature Molybdenum Element Furnace Control 
System" ] 


[Text] Introduction 


Among semiconductor manufacturing techniques, the sintering of raw ceramic 
needs to be done in a molybdenum element furnace at temperatures up to 1,750°C. 
Molybdenum element furnaces are characterized by high temperatures (up to 
1,800°C) and high power (tens of kw), and their control systems generally use 
the manually adjusted SCR output power method. Because the resistance of the 
molybdenum in the heating elements changes greatly as the temperature changes 
(as for example in the furnace we use when at 20°C the resistance is 0.0172, 
while at 1,750°C the resistance is 1.59), in order to control the temperature 
and safeguard the equipment, operations personnel must regularly adjust the 
output power. Within a furnace operation time of 12 hours, operations person- 
nel may not leave the premises. Nor can the accuracy of temperature control 
be guaranteed, as it oscillates +20°C. Because there is always a large flow 
of protecting H, gas in the furnace, this provides an unfortunate factor in 
controlling temperature. After the system in question used a microcomputer 
for control, there was a great improvement in performance, where single-point 
control accuracy reached 1,750°C +0.5°C, where the constant temperature area 
within 200 cm is +3°C, where the rate of product finishing has reached 99.5 
percent, and where at the same time labor intensity has been reduced. The 
system has excess temperature/cut-off thermocouple protection, and will sound 
an alarm. The rate of heat rise and temperature constant times may be set 
according to technical requirements, and the temperature can be maintained at 
any temperature. In addition, as conditions warrant, all set values may be 
changed without affecting program cperation. Temperature control is via the 
PID adjustment method, where to avoid integral saturation PI and PD disconnect- 
type control is used. Data collection is via digital filtering. To avoid the 
difficulty where changes in molybdenum element resistance is too great, the 
stepped change of mid-values method is used, which allows optimal control of 
the furnace throughout the temperature range. 











II. The System Hardware Composition 


System hardware is made up chiefly of a microcomputer, process input channel, 
process output channel, printer, digital clock timing circuit, alarm circuit, 
and power supply circuits. The microcomputer is a DJSO65A microcomputer, the 
CPU for which is the MBL6800, which has an instruction set compatible with the 
domestic DJS-060 series as well as with the imported 6800 series. Here, we 
have used only the single board computer part of the computer and the inter- 
face expansion board. The interface expansion board (I/OC board) has three 
parallel interface PIA MBL6821's, as well as the user parallel interface on 
the original board, which provides a total of four chips for the user to use. 
This computer has two MBM2716 EPROM's, for a total of 4K bytes; 2K bytes of 
the ROM is written with the monitor program, which leaves 2K bytes to be pro- 
grammed by the user. There are 4K bytes of static RAM altogether, of which 
1K bytes are used as stack for the monitor and the other 3K bytes are available 
for use. The system hardware flowchart is as shown in Figure 1. 


1. Process Input Channel 


The process input channel is outlined by the dotted line in Figure 1. This 
includes circuits for temperature sensing components, A/D converters, and 
isolating gates. 


a. This system is largely used at present for controlling a high temperature 
molybdenum element furnace, where the regular temperature is 1,750°C, the 
corresponding temperature sensor for which is tungsten-rhenium thermocoupler, 
where the thermocouple potential corresponding to temperatures of 0°C2#1,800°C 
is 02 28,030 microvolts. The change in thermocouple potential corresponding 

to each change of 1°C above 1,500°C is 13 microvolts. Because at the time 

this system was designed there was still no 16 bit A/D converters, we used the 
5 bit semi-digital voltmeter BS18A to substitute for an A/D converter. We used 
the 0.15 V part of the voltmeter, the least bit of which can detect 1 microvolt. 


b. Signal Gathering 


The output of the A/D converter is a 5 bit BCD coded signal, in all a 20-bit 
binary numeral, and when a parallel interface expansion board is added to the 
computer thena sample signal can be input once at the same time. By not using 
time shared gathering control circuits, but by adding a software digital filter, 
we have strengthened the reliability of the system. 


2. Process Output Channel 


As shown in Figure 1, this is composed of the PIA 2 B circuit, the PIA 3 A and 
B circuits, separator gates, D/A output and display, mid-value output and dis- 
play, point of inflection display, adder, power source and voltage compensation, 
and silicon controlled flip-flop circuits. 
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a. Mid-Value Output Channel 


Output in this system includes two parts: mid-value and D/A. In fact, mid- 
value is also D/A, and the D/A does not change for a particular segment of 
temperature area. Because the molybdenum element resistance will increase 100 
times from room temperature to 1,750°C, this is not like controlling an ordin- 
ary furnace where full power raises the temperature, but rather the method of 
depending upon changing mid-values is used to safely reach the temperature goal. 


It can be seen from Figure 1 that output is through channel B of PIA and that it 
is an 8-bit binary value, which means that it can divide the heating power into 
255 portions. Passing through the separator gate, one circuit drives an 8-bit 
LED, which displays the size of the output, and the other circuit connects to 
the mid-value D/A board. The D/A board is 12 bit, and according to need, we 
can connect the 8 bits of the PIA to the corresponding bits of the mid-value 
D/A. Because the 12-bit D/A is greater than the 8-bit output of the PIA, we 
can in this way flexibly select the size of the outpu.. After the mid-value 
output is added to the combination of the adder and the incoming D/A channel 
signal, the output then goes to the silicon controlled flip-flop circuit. Also, 
fluctuation of the power source voltage is added to the flip-flop circuit 
through hardware for use in high-speed compensation. 


b. D/A Output Channel 


Circuits in this portion are basically the same as those in the mid-value out- 
put channel, save for the fact that in the D/A output channel the 8-bit binary 
value output by circuit B of PIA 3 is the amount of adjustment output by the 
computer after the PID operation, ard normally that is fluctuating. Because 

the D/A board is 12 bit, power trimming can be as fine as 1/4095, which is bene- 
ficial to ensuring control accuracy. 


c. Point of Inflection Display Channel 


This outputs an 8-bit binary value through circuit A of PIA 3 and drives eight 
LED's through the separator gate. It displays which temperature areas are 
currently rising in temperature or are constant in temperature. This can dis- 
play the rising and stable temperature states at eight points of inflection. 


3. Interactive Equipment 
a. DJS-065A Microcomputer 


After the system is operational, beginning from the left of the computer display 
the first four positions will always display a hexadecimal value for the furnace 
temperature, which will correspond to the decimal value on the digital voltmeter. 
When the parameters need revising during operation, the STOP button may be 
pushed, which will temporarily halt operations. At this time, output will be 
maintained unchanging and the parameters can be revised or data can be examined. 
When finished, it begins again from the temporary halt program, and the system 
continues to operate. 








b. MX-80 III Printer 


We have selected the MX-80 model III dot matrix printer, and after the digital 
clock has generated a timing interrupt signal, the printer then prints out the 
millivolt value corresponding to temperature value at that time, and will at 
the same time print out the temperature curve. 


The timing signal can be set on the clutch plate from 1 to 60 minutes, and if 
arbitrary printing is needed the clutch plate can be moved to print out the 
millivolt values of the temperature for that moment. If printout is not needed, 
the printer switch can just be turned off. 


c. Other 


This system is also fitted with manual and automatic switching. When extra- 
ordinary conditions arise (such as excess temperature or thermocoupling), sound 
and light displays will automatically output that voltage has dropped to 

O volts. At this time, the situation can be handled by operations personnel by 
turning off the machine or switching to manual adjustment of output power. 


III. Mathematical Modeling of the System 


1. Adjusting Object Characteristics and a Flowchart for Controlling the System 
Transfer Function 


Ordinary electric heating furnaces are first-order inertial links, and this 
high temperature molybdenum element furnace is no exception, also being of the 
first order inertial link. Its transfer function is: 


W (CS) =K,e"''*/(1+T,S) (1) 


By measuring the rising curve of the furnace through experimentation, the 
amplification multiplier K, of the systemtransfer function can be found to be 
10, the time constant T, = 600 seconds, and the delay time ty; = 30 seconds. 


The system transfer function flowchart is shown in Figure 2. 
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Figure 2. Flowchart of System Transform Functions 
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In the figure, W.(S) is the input transfer function and is the group of broken 
lines determined by the process curve; W,(S) is the thermocouple change trans- 
fer function; Wp(S) is the power network voltage fluctuation compensation 
function; We(S) is the transfer function for controlling the target high 
temperature molybdenum element furnace; Wp (S) is the transfer function of the 
microcomputer's PID operation. 


The overall controlled output quantity is: 


Yey+y t+ ys (2) 


where y, is the mid-value quantity and a particular corresponding segment 
temperature does not change value. y; is the proportional feedback value 
determined by the hardware. The key here is to find the value y, in the PID 
operation of the microcomputer. 


2. Regulation Rules and Digital Filter for Sampling Data 


The ideal method of regulation for first-order inertial systems is proportion- 
integration-differentiation, that is, the PID regulation method, and the ideal 
calculation method by which to simulate regulation is: 


Y=K fe + (1/T,)fedt + T,(de/dt)] (3) 


The PID form of calculation for the discrete sampling equation is: 


Y.=Kyf{ eg #C1/T,) Beate, (Ae,/A7) | (4) 


Y, -- regulation value for PID K-order output 

Kp -- proportional coefficient 

T, -- integration time 

Tg -- differentiation time 

AT -- sampling time 

e, -- the deviation signal obtained by K-order sampling 
Ae,-e, 1 


where e;, is the deviation between the sampled measured value and the given 
value. In general, this value will be rather small when the system has entered 
operations in the vicinity of stable values. Therefore, interference cannot be 
neglected. Sampling is as in the following digital filter equation, which 
basically prohibits the effects of interference. We proceed this way with one 
signal; thatis, we successively sample five times, then eliminate the largest 
and the smallest, then take the average value of the three remaining mid-values 
as the value at this time. But the difference e, between this value and the 
set value is not taken as the deviation value at this sampling. Rather, the 
deviations of the preceding three occasions are averaged with this deviation 
and substituted for the deviation value at this time. Showing ey, we get 


= (e, +e +e )/4 (5) 


e, mwa * Sue * Su3 
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de, /AT = (e, - e, , + 3e, _, - 3e, _,)/6AT (6) 


/ST is the average value for deviation differentiation, so substituting 
first (5) then (6) in equation (4), we get: 


Y,= K,{ Ce. Heys Hey-g HOn_3)/4 


7 7 
+(4T/T1) Sey + Tele, eyes ' 


+3ey-, -3e,-1)/6AT} 


where Kp» T,, and Tg can be roughly set up beforehand through Kr, Te, and Tr of 
the furnace transfer function as follows: 


T, =2r, =60; T, =0.6r, = 15; 
K, =T,/(0.85K,t,) =2.15 


Accurate adjustments require that in actual practice optimal parameters be 
obtained, so at 1,750°C our final adjusted values are: 


x, = 4; TT, = 80 seconds; Ty = 20 seconds, 


IV. Applications Software 
1. Main Temperature Control Routines 


The flowchart for the main routines for controlling temperature is shown in 
Figure 3, which we now explain in the following points. 


a. After PIA initialization and clearing of the working elements to 0, there 

is a delay of 2 seconds to allow the entire regulation period to be 3.2 seconds. 
This lets it match the control target and will also leave some time when in the 
future we control several furnaces. 


b. Sampling of the temperature millivolt values is continued for five times 
(separation being 200 ms because the sampling period of the A/D converter is 
200 ms), and after these five values have been digitally filtered they are used 
for identification. If this value is very small or equal to 0 it indicates that 
the thermocouple has already burned off or that there has been some other mal- 
function. The mid-values will be turned off, the D/A converter will be turned 
off, and an alarm will be sounded. If this value is greater than the number of 
millivolts in the thermocouple corresponding to 1,800°C, then this indicates 
excess temperature, and the mid-values and D/A will be turned off immediately, 
and an alarm sounded. 


c. The final routines determine in order whether 1,500°C?; 900°C?; 600°C?; 
or 300°C? If the temperature is lower than 300°C, the production process does 
not need a rising temperature slope below 30°C, and also because at this time 
the furnace molybdenum resistance is quite small. To reliably raise the 
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temperature, the method of adding $03 to the D/A every 20 seconds is used to 
gradually increase the output power. When the temperature has exceeded 300°C, 
the entire temperature process is divided into four parts, which allows each 
segment to correspond to different mid-values (stored as element $B) and D/A 
(stored as element $A), ensuring that the power will not become too great and 
burn components. 


d. Determines to which segment a temperature corresponds, then after filling 
in the parameters for this segment in the working element, by distinguishing 
whether the high bit of control register BCOl is 0 or 1, it is determined 
whether it is necessary to raise the temperature or to maintain it. If an 0, 
it indicates that there is no pulse input to jump to the raise temperature 
routine. If al, it indicates that there is a puise input and processing is 
to be as needed for maintaining the temperature. Please note that at the end 
of the routines to raise temperature, they return to the temperature mainte- 
nance routines. 


e. After jumping to the temperature maintenance routines, it is first deter- 
mined whether the current temperature is lower than the low limit temperature 
S10w that has been set. If it is lower than 9),,,, then the mid-values and D/A 
are turned completely on (the mid-values and the D/A of the corresponding 
segment are not the same), which allows a rise in temperature at full speed; 
when the current temperature is higher than the upper limit temperature Sigh 
that has been set, then the mid-values and D/A are turned completely off for 

a reduction of temperature at the greatest speed. 


Only when the temperature is between the set values Ojo9y and Onigh can the 
routines enter the normal range for controlling temperature, solving for the 
deviation value ex. If the absolute value of ey is greater than 5°C, PD con- 
trol is used; if less than 5°C, PI control is used, which allows for a better 
control process. The control value Yy, is derived from PID calculations. It 
is used to regulate the furnace heater flow I[,, the total heating current of 
the furnace being composed of three parts: 


T= 1, + Ilo + Is (8) 


where I, is the mid-value current, corresponding to the temperature within a 
particular range, this current does =ot change. I») is the compensation cur- 
rent as determinedby the hardware and the power source voltage reciprocal. 

I, is the control current as calculated by the sicrocomputer through PID cal- 
culations. Cwrresponding to any particular temperature, when the deviation 
e, is 0, any value may be present. Above 1,500°C, the system routines are 
preset to 88, and when this value is added to or subtracted from the calcu- 
lated Yue the result is sent to the D/A, which is the I, that we have derived. 


f. The raise temperature routines rise linearly according to the set value, 
the furnace temperature reached by the principle of rising. The purpose is to 
allow the rising temperature curve to fit the process requirements. 
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2. Printer Routines 


This temperature control system is provided with an MX-80 III printer, a dot 
matrix printer that is a product of the Japanese Epson Co. and that has been 
modified by the Tianjin Hongxing Plant. Printing is by means of interrupts 
generated to the computer. This computer provides the user with a screen 
interrupt IRQ (active low), as shown in Figure 1, and the signal to print is 
a timing signal (adjustable from 1 to 60 minutes) generated by the digital 
clock. Output generates an interrupt from the PIA 4 CAl. 
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Figure 4. Flowchart of Print Routines 
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A flowchart for the printing routine: is shown in Figure 4, some points about 
which we explain here: 


a. The printercan only receive ASCII code, so after the binary values have 
been changed to decimal, they must still be changed to ASCII before outputting 
to the printer. 


b. Because we are using compressed printing, there is a maximum of 132 char- 
acters [in a line], while the millivolt values for the temperature are five 
place decimal values, which leaves 127 character positions. We have allowed 
the largest available space to be 125, the hexadecimal value »f which corre- 
sponding routine is $99C4. 


c. The printing format is to first print out the decimal value of the milli- 
volt value of the current temperature, after which is printed an asterisk. 
The position of the asterisk is determined by the following formula: 


M= (T - A)/C (9) 


where M -- the output space value 
T -- printed temperature microvolt value 
A -- initial microvolt value 
C -- the microvolt value represented by each space. These routines 
assume 20 microvolts. 


We can see from equation (9) that there will be an asterisk printed after M 
number of spaces. 


The initial value is chosen by the operations personnel, where it is placed at 
addresses $0210 and $0211, and if less than the initial value or exceeds 125 
spaces, then it is not printed. 


The printing routines inventory has been omitted. 

V. Conclusions 

The temperature controlling accuracy of this terminal is +0.5°C (600¥1,750°C); 
it can run for long periods. With slight modifications, this software and 
hardware can control several devices by itself. But as microcomputer prices 


get lower and lower, from the point of view of highest reliability we still 
hope to develop microcomputer single installation control systems. 
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APPLIED SCIENCES 


MICROCOMPUTER CONTROLLED ANALOG TESTING SYSTEM DESCRIBED 


Beijing DIANZI JISHU YINGYONG [APPLICATION OF ELECTRONIC TECHNOLOGY) in 
Chinese Wo 2, 25 Feb 86 pp 18-20 


[Article ty Wei Fengwei (7614 7685 5633]: "Microcomputer-Controlled Analog 
Quantity Testing System") 


[Text] In any special field there mast be accurate, error-free testing of 
analog electric quantities, and this must be quick and the instrumentation 
light and convenient. This is quite a high requirement for automated test 
equipment. By using microcomputers in the testing and control of analog 
quantities we can easily realize this requirement. Regarding the 
microcomputers, all they need are some hardware interfaces together with the 
corresponding testing software and they can then quickly and accurately 
accomplish the testing requirement in accordance with pre-arranged routines. 
They can also generate an alarm signal when volt or current have exceeded a 
standard. They can provide the location of the error, and the test results 
can be completely output to display or printer for file storage and analysis. 


I. A Brief Description of the Testing System 
The testing system is composed of the five parts that are the object to be 


tested, the controller instrument that is microcomputer based, the test 
instruments, a power supply, and output equipment. A diagram of the analog 


quantity test system is shown in figure |. 
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Figure 1 Flowchart of Analog Quantities Test System 








Control commands are generated by the user program in the EPROM or the 
microcomputer, which allows automated testing of multiple circuits and various 
Kinds of signals. Each signal that is tested comes in through cables to the 
test instrument's signal normalized circuit, gated through multichannel analog 
switches and gates, sequentially sent into ac/dc transformer circuits and dual 
integrated A/D converters, and is then sent to the microcomputer via an input 
interface where it undergoes data processing before being output for display, 
when the test results are printed out. A partial power source schematic for 
this automated test system is given in figure 2. 
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Figure 2 Power Supply Schematic 


a. The obdject of testing 


The objects of testing can be electric signals that are the output from power 
supplies for electronic equipment and from voltage, current, and signal 
sources, as well as working points on an electronic circuit; all output 
current not less than 50 microamps; DC voltage signals of 0+ or - 40 volts; 
effective values of 25 volts of low frequency sinusoidal voltage; peak values 
of 40 volts of low frequency square waves, all of which can be tested quite 








accurately in this system. Length of the transmission cable is determined by 
actual conditions, the only requirement being that the internal resistance be 
as low as possible. Theoretically, the voltage values to be tested may be 
raised without affecting the tested accuracy. But if the tested source load 
capacity is quite weak, during the process of transmission the signal will be 
affected by interference, which will affect the accuracy of the test results. 


b. The controller 


The controller is composed of three parts: the microcomputer, a sampling 
decoder switch matrix, and an input interface. 


The microcomputer used is a Z-8000 16-bit single board computer. But the TP- 
801 8-bit computer can also be used. 


Sampling control is achieved through the 8-bit data lines provided by the 
microcomputer. The D6 and D7 among them are the status control bits, which 
two bits can decode four different codes. It can select four different 
working states of AC-DC transformer circuits: afterpulse, maintain peak value, 
Maintain negative peak value, and maintain. Lines DO through D5 are the 
sampling gate controller bits, through which can be decoded 64 different 
codes, used for controlling an 8 X 8 switch square matrix. There can be 
automatic gating and testing of signals on 64 circuits. The testing system 
described in this paper actually uses 40 circuits, the other 2. being for 
future use. The switch matrix uses the non-contacting electronic multiplexing 
analog switch CD4051. 


The input interface is a 16-bit parallel input interface. Two levels of 
transfer have been designed for this system, with four 4D flip-flop 74L8175's 
acting as data buffers; two 8-bit tri-state latches, 74L8244, latch the 
computer's data bus. 


c. The test instrument 


The test instrument is made up of a normalized circuit, AC-DC transformer 
circuits, and an A/D converter. 


The normalized circuit is actually a system input circuit that has been 
designed as an RC network with filter, current-voltage converter, and voltage 
level conversion functions. By adjusting the R1-2 ratio we can select the 
sampling ratio, allowing the input signal to uniformly become a voltage signal 
within a particular amplitude, to prepare for gating and testing. 


The A/D converter is the key link in the test instrument. We have used the 
dual integrated A/D converter, 5G14433. This device is a product of the 
Shanghai Component Factory #5, and outputs BCD code. It uses analog zeroing, 
digital zeroing, repeated analog zeroing, integration of input voltage, and 
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accomplishes first-order conversion of analog quantities-digital quantities 
for reference voltage integration, integrator constant-current discharge, and 
for counting within the current discharge time. R2 in figure 2 is used to 
adjust components connected to the clock frequency of the converter, as_ the 
clock frequency determines the speed of the A/D conversion. The conversion 
speed of this system is 4 times per second. 


The AC-DC converter is a dual integrated A/D converter 5G14433 designed 
advance circuit. It function is to convert the normalized voltage signal to 
be a DC voltage signal of less than 2 volts, to prepare for modular 
conversion. The working process for the AC-DC converter is controlled by the 
computer, and four working states are automatically selected according to the 
requirements of the test signal. 

d. Power supply 


For the power supply a floating, dual shielded secondary power source is used. 
e. Output equipment 

1. A digital display device on the microcomputer; 2. a printer. 

f. Functional set-up of the system 


1. Sampling of the input signal is gated through control of the 
microcomputer, automatically gated according to the user-designed program. 


2. Working states oof the AC-DC converters are controlled by the 
microcomputer, automatically programmed to switch over. 


3. Various reset coefficients and standard values for tested electrical 
quantities are programmed into the EPROM, so data handling proceeds 
automatically. 


4, Actual values of tested signals can be read directly from the display 
device or from the printed results, the dimensions being identical with 
reality. 

5. Automatic adjustment of zeroing. 

6. Automatic initialization and both manual and initial reset. 

Il. A Flowchart of the Test Routines 


A flowchart of the test routines is given in figure 3. 
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Figure 3 Flowchart of Testing 


III. Further Developments of System Applications 


What we have discussed above is an automatic testing system controlled by a 
microcomputer. It can also be used for multimonitoring or relative constant 
control of physical quantities such as temperature, speed, and water pressure. 
As long as each monitoring point (or control point) of the object to be tested 
is equipped with sampling devices and electrical quantity converting devices 
that can convert the physical quantities to electrical quantities, it will 
work. The flowchart of the routines used for monitoring is shown in figure 4. 
The portion indicated by dotted lines in the figure are the places in the 
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program process and monitor testing program process that are not the same when 
used for relative constant control, the identical portion having been 
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Figure 4 Flowchart of Monitoring Routines 
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IV. System Anti-Interference Measures 


Whether it is the testing of analog quantities or the process of monitoring or 
controlling, when examining the strengths and weaknesses of a system anti- 
interference capability is an important index. If the system is operating 
under harsh conditions, the relative strength of its anti-interference 
capability not only directly affects the reliability and accuracy of output 
results, but it even concerns the success of the system. To improve anti- 
interference capabilities in this system, we have taken the following 
measures: 


a. We have chosen a dual integrated A/D converter that is highly resistant to 
interference; 


b. Filtering links have been added to input circuits to prevent random high 
frequency pulse interference; 


C. The microcomputer and test instruments both have high level shielded, 
floating secondary power supplies and automatically loop back to prevent 
common mode interference; 


d. The differential input method is used to prevent ordinary state 
interference; 


e. The transmission cable has twisted-pair shielded wire, and the test 
circuits are shielded, all of which prevents cross interference. 


Experiments have shown that anti-interference measures used in this system are 
effective. When transmitting over a 150 meter line, test accuracy of the 
system is not degraded. 


12586 
CSO: 4008/1075 
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[Text] I. Introduction 


During the operation of a nuclear power plant, most of the released 
radioactive elements with the exception of the low-poison group such as 
inert gases and tritium, are contained in th radioactive solid wastes. 
Therefore, management of radioactive solid waste is an issue of great 
importance. 


The treatment of radioactive waste involves converting the waste materials 
into a form that is suitable for shipping, storage and disposal (including 
condensation, compression, solidification and packaging). The disposal of 
radioactive waste is to properly seal and store the waste so that before 
the radioactive elements decay to a sufficiently low level, they will be 
isolated from the human and biological cycles, and will never cause harm to 
human beings. 


In 1983, the International Atomic Organization held a conference on 
radioactive waste management in Seattle, where 31 countries and 9 intarnational 
organizations participated. It was generally believed that the basic technical 
problems of disposing low-radiation level waste materials had been solved, 

but puyeaee improvements are required to develop a complete waste disposal 
system’. 


II. Unique Features of Radioactive Solid Waste from Nuclear Power Plant 


It should be pointed out that the radioactive solid waste from nuclear power 
plants does not include the exhausted fuel elements from the reactor or the 
high-radiation glass byproduct from processing exhausted fuel elements. A 
small amount of high-radiation waste such as replaced reactor parts are stored 
with the exhausted fuel elements where they undergo decay and then shipped 

out of the plant. 
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A small amount of solid waste (e.g., contaminated equipment and parts) 

may be treated as non-radioactive waste or reused after decontamination or 
storage and decay. Most of the solid wastes (e.g., solidified residual 
liquid from evaporation, filtered slush, filter core, wasted resin, burned 
ash, or other types of dry waste) contain more than 90 percent of the 
radioactive elements. These solid wastes are primarily low to medium 
radioactive elements with relatively short half lives and exhibiting 8 - y 
radiation; the content of super uranium, long half-life elements exhibiting 
a radiation is very low. The level of activated corrosive products sygh 
as 58Co (Ts = 71 days), ®%Co (Ts = 5.3 years), 2>Fe (Ts = 2.7 years), Fe 
(Ts = 45 days) is quite high, but because of their short half lives, generally 
they will have no adverse effect after storage for 50-100 years. The 
radioactivity after storage primarily comes from fission products 99Sr (Ts = 
28.5 years) and 137cs (Ts = 30.2 years). 


The annual output of solid waste from a nuclear power plant depends on the 
reactor power, the type of reactors, the operation and management of the 
power plant, the condition of maintenance and repair, and the waste treatment 
techniques (27, The amount of solid radwastes produced by a 1 million kW 
light-water reactor #* the U.S. is 400-1300 m3 per year [ 1; a similar 
reactor in Japan produces (600-1000 m? solid radwastes per year [4]; and 

a 900,000 kW PWR of the French Framatome Atomic Energy Co. produces 400-500 m3 
solid radwastes per year. 


The above discussion shows that the unique features of solid radwastes 

from nuclear power plants are: large volume (a 1l-million kW PWR power 
plant produces 400-600 m? solid radwastes per year; a BWR power plant 

with the same reactor power will produce 70 percent more), and low specific 
radioactivity, most of which is attributed to short half-life nuclear 
elements with 8 - y radiation. According to the law of radioactive decay, 
the level of radioactivity is reduced by a factor of 1,000 after 10 half 
lives and by a factor of 1 million after 20 half lives. In other words, the 
radioactive waste must be isolated from human and bilogical cycles for 
300-600 years. 


III. Ocean Disposal 
Ocean disposal of solid radwastes has the following merits: 
1) The ocean is isolated from human en. ‘ronment; 


2) The structure of the ocean floor used for waste disposal is geologically 
stable; 


3) Because of the lack of free oxygen near the ocean floor, the waste 
containers can remain intact for a long period of time without suffering 
oxidation; 


4) The ocean emits a high level of natural radiation, even if there are leaks 
of radioactive elements from the disposed wastes, they will be absorbed in 
the fluctuations of natural radiation. 
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The Nuclear Energy Office of the Organization of Economic Cooperation and 
Development had studied the problem of ocean disposal of low radiation 
solid wastes, and had established guidelines and clearly defined monitoring 
responsibilities. Disposals of nuclear wastes in the Atlantic Ocean had 
been monitored by an international organization; specifically, at a depth of 
5,000 m, five member nations had disposed 35,000 waste containers with a 
total weight of 11,000 tons and total radiation level of 8,000 curies. 
Subsequently, eight European countries (France, West Germany, Belgium, 
Italy, Holland, Switzerland, Sweden and Great Britain) participated in 11 
ocean disposals; the average depth of disposal was 4,500 m [5]. 


In order to prevent the radioactive wastes from polluting the ocean, a 
conference was held in 1972 in London to discuss a treaty on ocean disposal. 
The treaty requires that each country establish a national-level organization 
to control the ocean disposal of wastes, ban the disposal of high-radiation 
wastes, and that disposal of low and medium radiation wastes must be approved 
and managed by a special office. It also establishes specific regulations 
for selecting the disposal sites, and requires that they must be at least 
4,000 m deep, and must away from continental mounts, deep-sea resources 
and underwater cables 


Based on an assessment of the radioactive ecology of disposed radwastes in 
deep ocean regions, the International Atomic Energy Organization has banned 
the ocean disposal of high-radiation wastes; all radioactive elements are 
classified, and those which exceed the threshold are banned from occan 
disposal. An inspection of the disposal site in the northeast Atlantic 
Ocean shows that it meets the requirements established by the International 
Atomic Energy Organization. 


Between 1946 and 1970, the United States opened more than 50 ocean disposal 
sites, 90 percent of the waste materials and 95 percent of the radioactive 
elements are concentrated in four of these sites, but they did not meet the 
specifications of the London treaty or the Organization of Economic 
Cooperation and Development. Since the method of shallow ground burial 

of solid radwastes was successfully developed by the Department of Energy, 
in 1970 the United States decided to stop ocean disposal of radioactive 
wastes in order to protect the ocean environment. In an environmental 
assessment of the four disposal sites, containers were retrieved from the 
ocean floor ranging from 2,800 m to 3,800 m in depth; careful inspection 
showed that the containers were compressed toward the center, and there 
were small amounts of radioactive elements surrounding the container; however, 
there was no evidence of damage to human beings or the ocean environment [7], 


In 1949, Great Britain began disposing packaged containers of low- 
radiation solid wastes in the Atlantic Ocean; between 1967 and 1971, it 
participated in the experimental joint ocean disposal project sponsored by 
the Nuclear energy Office of the Organization of Economic Cooperation and 
Development, and accepted international monitoring. By the end of 1983, 
the total amount of ocean disposed waste materials reached 77.53 tons. 

They believe that ocean disposal will not lead to detectable danger, and 
therefore will continue to do so. They plan to invest in the renovation 

of disposal ships and installation of satellite navigation equipment; the 
amount of disposed radwaste is expected to increase at a rate of 20 percent 


per year. 
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In order to assess the safety of ocean disposal, J.upan conducted high-pressure 
filtration tests and container integrity tests undi:r simultated 5,000-m deep 
ocean conditions. It also conducted safety inspections of disposal sites. 


In 1983, the London Waste Disposal Treaty Conference passed a resolution to 
ban ocean disposal of radioactive wastes, but this resolution had no real 
restraining authority. Protecting the ocean environment and preventing 
radioactive pollution has become a problem that is of concern to every 
nation in the world; a final resolution is expected tw be reached during 

an international conference in 1985. 


IV. Land Disposal 


Land disposal of radioactive —_ ree is an effective method to protect 
human beings and the environment [8 ; it accounted for 99 percent of the 
radwaste disposals in the past (9), ‘There are three different approaches to 
land disposal. 


1. Shallow Ground Burial 


The shallow ground burial approach is to bury the radioactive solid wastes 

in a trench which is 30-240 m long, 6-12 m wide and 6 m deep. The waste 
container is first covered with 1-3m of dirt then with a layer of clay; a 
stone or metal market will be erected to indicate the volume and radiation 
intensity of the disposed wastes. The United States has 6 commercial disposal 
sites which can accommodate 764,000 m3 of solid wastes. Currently, all the 
low-radiation solid wastes produced by the nuclear power plants in the U.S. 
are disposed using the shallow ground burial approach. At some of the earlier 
burial sites, because of the crude construction and site selection, and poor 
packaging and management, there was evidence of leakage of radioactive 
elements. Upon detection by the monitoring system, the disposal sites were 
closed down immediately to avoid further damaging effects on the public and 
on the environment. 


The construction and utilization of shallow disposal sites can be divided into 
four stages: 2-3 years of building the site facilities and digging the trenches; 
15 years of disposal service; 150 years of monitoring and limiting land use 

after closing down the disposal site; after that the land can be used without 

any constraint. Because of the different conditions of hydrology and geology, 
the cost of building a disposal site differs considerably from one country to 
another. 


In Great Britain, the Nuclear Industry Radioactive Waste Commission divides 
the waste materials designated for shallow ground burial into different 
categories; low-radiation wastes can be directly placed in storage tanks; 
medium-radiation waste tanks must be shielded with cement before shipping; 
waste materials with high 8-y radiation are stored in tanks and shipped to 
the disposal sites in shielded conta‘ners. The design guidelines for the 
burial trenches are: 


a) selecting a trench site with the proper geological structure of with clay 


soil; 
b) the sidewalls and bottom of the trench must be inclined; 
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c) water drainage ditches must be constructed around the burial trench; 
d) monitoring holes must be installed in the vicinity of the trenches. 


The first type of burial trench is 150 m long, 8 m wide and 6-9 m deep. the 
bottom of the trench is inclined toward the water pit, and is covered with a 
layer of l-m thick gravel and sand for water drainage. When the waste 
materials reach a height of 3-5 m, they are covered with clay to minimize 
water infiltration. 


The second or third type of trench is 70 m long, 20 m wide and 19 m deep. 

The bottom of the trench is inclined toward the water pit, and covered with a 
layer of porous cement; the water pit is equipped with monitoring tubes, and 
the trench walls are constructed with steel reinforced concrete. The shielded 
waste tanks are lowered into the trench using a fork lift; after reaching 

a height of 8 m, they are covered with a l-m thick steel-reinforced concrete 
and coated with a layer of water-resistant asphalt, then filled with a layer 

of 3 m-thick packed clay. On top of the trench is a 1 m-thick steel 
reinforced concrete structure covered with a layer of of 5 m-thick packed clay 
to prevent rain water from entering. 


In France, the disposal of solid radwastes from nuclear power plants is the 
responsibility of a national radwaste management bureau. The [wangshen] 
radwaste storage pool located in northwest France has a storage capacity of 
400,000 m?; since storage began in 1969, 200,000 m3 of it has been filled. 
In France, the solid radwastes produced by a nuclear power plant are stored 
inside the plant for 15 days before being shipped to the national waste pool 
for disposal. 


2. Storage in Abandoned Mines (or Caves) 


This form of storage is primarily used by small, densely populated European 
countries with stable geological structure. The abandoned [Hesai] salt mine 
in the northeast part of West Germany is free of ground water and geologically 
stable; it has an area of 4.2 square kilometers. In a region 490-750 m below 
the surface, 130 chambers were excavated with a total capacity of 3.6x10° m>; 
they can accommodate 140,000 barrels of low-radiation wastes and 1,300 

barrels of medium-radiation wastes. 


Abandoned mines generally have sufficient depth to have any effect on human 
activities. But they were designed and constructed to meet mining 
requirements; whether they are safe for disposal of radioactive wastes remains 
to be investigated. 


3. Storage at the Plant Site 


In technical report No. 198 of the International Atomic Energy Organization, 
"Guidelines for Safe Treatment of Radioactive Wastes of Nuclear Power Plant", 

the feasibility of plant site storage was confirmed. The guidelines stated 

that "for a reactor location with favorable hydrological and geological 
conditions, it is feasible to use shallow ground burial to dispose certain 
radwastes at the plant site." Japanhas accumulated 460,000 barrels of solid 
radwastes stored at the reactor sites. The United States had also suggested 
plant site disposal plans because of increasing disposal costs due to 

shut-down of three of the four disposal sites in the eastern part of the country. 
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Final disposal of the solid radwastes and the plant will take place when 
the nuclear power plant reaches retirement age. Great Britain on the other 
hand, believes that plant site storage of solid wastes can be considered a 
form of final disposal. 


The long-term process of radwaste disposal may be affected by natural 

processes or by human activities, which cause radioactive elements to enter 

the human environment through water infiltration and diffusion. Because of the 
difficulty in finding an ideal disposal site, man-made screens (e.g., engineering 
structures) are often used to enhance the isolation effect and the depth of 
burial is increased to improve safety. The disposed solid wastes must meet 

the following requirements: 


a) there must not be any leakage of liquid or gas; the U.S. Nuclear 
Regulatory Commission specifies that the amount of free liquid in a waste 
tank must be less than 1 percent; 


b) the container must have sufficient strength to withstand loading, 
unloading and docking operations; generally, 200-liter steel tanks with wall 
thickness of at least 1.5 mm are used; Japan specifies that the waste tanks 
mist be able to withstand a load of 5 tons; the United States specifies that 
the containers must remain corrosion resistant for 20 years and they must 
be retrievable; 





c) the content of super uranium a radiation in the disposed wastes must be 
strictly controlled. 


Solid radwastes of nuclear power plants have much lower screening 

requirements for the disposal sites than high-radiation, long half-life 

waste materials. They can be disposed in shallow or medium-depth trenches, and 
that the heat generated due to decay can be ignored during the disposal period. 


V. Suggestions 


Disposal of solid radioactive waste of nuclear power plant is an important 
part of waste management, and directly affects the cost of waste management. 
During the disposal period, the effect of radiation on the environment and on 
human beings mist be kept "as low as feasible." Disposal of solid radwaste 
not only is a costly task but also requires the allocation of a large amount 
of land; therefore, reducing the volume of the waste material is a critical 
issue of waste management in countries with high developed nuclear power 
capability as well as in developing countries. 


Disposal of solid radioactive waste of nuclear power plant is a complex 
problem which has technical, economic and social implications. In view of 
the unique characteristics and far-reaching effects of radioactive wastes, 

any policy concerning the disposal of radwastes must be proclaimed in the form 
of national laws or regulations so they will be enforced with sufficient 
authority. A national radioactive waste management organization should be set 
up to establish the necessary regulations and standards including radwaste 
management procedures, waste classification standards, and performance test 
standards; it should also be responsible for such tasks as mapping the 
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hydrological and geological structures of the earth stratum, studying the 
mechanism of migration of nuclear elements, selecting appropriate disposal 
sites for land burial in order to support China's nuclear power industry. 
Environmental protection is a basic policy of this country; therefore, economic 
development and urban/rural development must be coordinated with environmental 
de ‘lopment in terms of planning and implementation in order to reach a 
bal..uce between economic benefits, social benefits and environmental benefits. 


3012/9869 
CSO: 4008/77 
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[Text] I. Introduction 


The U0> elements in a light water reactor (U02-LWR) produce approximately 

1 percent industrial plutonium after irradiation; most of the plutonium 
produced are 239Py and 24lpu, which can be used as fuel for theremal neutron 
fission. Since the 1950's, western countries in Europe and North America 
have conducted studies and industrial tests of the recycled utilization of 
treated nuclear residues of UO?-IWR. In 1957, the first manufacturing plant 
of mixed uranium and plutonium oxides (MOX) began operation; in 1963, the 
first batch of experimental MOX elements were tested with BR, irradiation 

in a reactor. By January 1984, 13 PWR's and 12 BWR's were using MOX elements 
with a maximum burnout value of 71 Gwo/tul!!. 


Since the beginning of the LWR industry, recycling U and Pu has been a topic 
of great interest because of the low utilization rate of LWR fuel and the 
limited uranium resources. While building FBR is one way to solve this problem, 
a completely new industrial system which has more stringent requirement 

than that of LWR mst be developed. It is estimated that such a system will 
not be commercially available until around 2025, and even then its initial 
economic benefits will not be able to compete with LWR's. Consequently, 
suggestions have been made to use a MOX-LWR fuel cycle strategy during the 
LWR and LWR-FBR transition period, i.e., using the retrieved uranium and 
plutonium from LWR post-processing to make MOX-LWR elements, which can then be 
reused in LWR's. The author believes this to be an effective approach to 
improve the fuel utilization LWR, to develop new sources of fuel, and to 

speed up the accumulation of plutonium for FBR's. 


II. Recycling of Uranium and Plutonium in LWR's 


1. Recycling of Plutonium 


The industrial-grade plutonium (also called primary plutonium) produced from 
burning UO>-LWR in LWR's amount to approximately 0.9-1 percent of the total 
reserves of heavy metals. [It contains approximately 1 percent of 738Pu, 58-68 
percent of 239pu, 20-30 percent of 240py, 9-10 percent of 24\pu, and 2-6 percent 
of 242py, Primary plutonium can be mixed with depleted uranium, natural 

uranium and retrieved uranium from post-processing to produce MOX 
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elements, which are reusable in IWR's. Energy will be released from thermal 
neutron fission of 239py, 241py and 235y, and the neutrons will again be 
mega by 238y, 238Ppu and 240py to produce fissionable plutonium (239Pu, 
24lpy). Therefore, by burning recycled eclements, there is no intrinsic or 
quantitative reduction in the fissionable plutonium; in other words, there is 
no intrinsic change in fissionable contents between secondary plutonium (from 
MOX elements) and primary plutonium, or between tertiary plutonium and 
secondary plutonium. This is the basic difference between plutonium and 

235y in the U0, fuel elements. 





Fig. 1 shows the composition of plutonium isotopes at various stages of 
recycling. After the fifth recycle, the proportions of plutonium isotopes 
remain unchanged; they do not vary with the number of recycles or with the 
type of uranium used (natural uranium, depleted uranium, or retrieved 
uranium from post-processing). The “equilibrium” isotope contents are 
shown in Table 112), 


Fig. 2 shows a plot of the variation of fissionable plutonium (Pufje) as 4 
percentage of the total plutonium with the number of refuelings in the 
reactor. Initially, Pugig is about 77 percent of the total plutonium, then 
it gradually decreases and stabilizes at about 60 percent after 20 

refuelings. This “equilibrium composition” of industrial plutonium is totally 
acceptable for maintaining the production of MOX-LWR elements. 


Table 1. Pu Isotope Content After Fifth Recycle of MOX Elements 











Pu isotopes 7% 
- §___4_,_3_,_#_, 


No. of refuelings a i. 





No. of refuelings 


Fig. 1. Composition of Recycled Pu Isotopes in MOX Flements 
Key: 
(1) IR: lst Pu recycle 
(2) SR: Sth Pu recycle 
(3) equilibrium composition of Pu isotopes, percent 
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Fig. 2. Variation of the Proportion of Fissionable Plutonium with Number 
of Refuelings 


Key: 
1. use of natural uranium as basic fuel 
2. weight of Pufis in MOX elements 


The key considerations for optimum recycling of industrial plutonium are 

the following: a) the plutonium should be allowed as much time as possible 
in the reactor to maximize energy production; b) the elements should have 

an adequate cooling period after being removed from the reactor to ensure 
disintegration of the short-half-file fission elements in order to reduce the 
level of irradiation during post-processing; c) the 241lPu should be utilized 
as much as possible in the reactor before significant decay (T1=14.4 years) 
occurs. 





Based on the above considerations, the Nuclear Industry Service Co. (NIS) of 
West Germany suggested that the safest and most economical approach is to 
have a 2-year period outside the reactor, which includes 1 year of cooling 
and decay and 1 year of post-processing and element production. 


The standard operation for both PWR's and BWR's is to have one refueling 
per year; in PWR, one-third of the fuel is replaced each time, and in BWR, 
one-fourth of the fuel is replaced, hence the plutonium recycling periods 
are respectively 5 years and 6 years. If the life of a reactor is 30 years, 
then PWR's can be recycled 6 times, BWR's can be recycled 5 times. 


The initial proportion of MOX elements in the reactor core is only 3 percent 
for BWR and 4 percent for PWR. As the number of cycles increases, plutonium 
production also increases; by the 20th refueling, plutonium consumption and 
production reach an equilibrium, i.e., the so-called SGR (self-generated 
recycling) condition, and the proportion of MOX elements becomes approximately 
30 percent. Table 2 shows the data of fuel elements for different 

refuelings; the data do not vary significantly if natural uranium is replaced 
by depleted uranium or post-processed uranium, 


Table 2. Fuel Element Data in a Reactor (Natural Uranium) 
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Based on the “optimum self-generated recycling; mode of operation, plutonium 
remains inside the reactor for 20 or 18 years of the 30-year reactor life, i.e., 
more than three-fifths of the time it is producing energy for power generation. 
However, if the current “delayed plutonium recycling" mode of operation is 
used, during each cycle it remains outside the reactor for 7 years, plus 1 
additional year required for post-processing and element production, i.e., the 
period of storage and processing outside the reactor is four-fifths of the 
reactor life. Clearly this is undesirable from the point of view of both 
economy and safety. 


2. Recycling of Uranium 


Most of the uranium in the U0? elements from a LWR contain 238y and 0.8-0.9 
percent of 235y, both of which can be returned to the LWR without much 
difficulty, but there are also accumulations of 232y, 234y and 236 in the 
retrieved uranium. Although the 232U content is very low, its strong y 
radiation requires special protective measures and the necessary modifications 
of the industrial flow process [3], 234y and 236U have strong neutron- 
absorption capability, hence the presence of 234y and 236 tend to reduce the 
reactivity of the elements, whick can only be compensated by increasing the 
fuel content. 


In terms of 234y and 236, the first uranium recycle can save 15 percent of 
natural uranium, but the second recycle can only save 3 percent of natural 
uranium, hence its economic benefits are limited. 


The above discussion shows that from the point of view of cost effectiveness 
and minimizing processing difficulties, uranium should be recycled only once. 


3. Experience of Using MOX Elements in LWR 
There are generally two types of MOX modules used in a LWR: 


a. "Plutonium Island" Type Modules. In this case, the MOX elements are 
placed in the center of the reactor core, with enriched U0) elements placed 
around the periphery. This type of module is suitable for BWR's or reactors 
with large water gaps between the modules. 


b. "All Plutonium" Type Modules. This type of module contains only MOX 
fuel elements and is used primarily in PWR's. In this type of reactor, the 
effects of guide tubes on flux and peak value coefficients can be corrected 
by properly selecting the enrichness distribution diagram. 


In a typical PWR module, approximately 60 percent of the MOX elements are 
highly enriched (3.9-4.6 percent fissionable materials), 8 percent are slightly 
enriched (2.4-2.8 percent fissionable materials); and the remainder are 
moderately enriched (3.2-3.9 percent fissionable materials). 


In a typical BWR module, approximately 50 percent of the MOX elements are 
highly enriched (3.2-4.6 percent fissionable materials), the other 50 percent 
are slightly enriched (2.6-3 percent fissionable materials). 














With regard to reactor core management, one can follow a procedure similar 
to that used in U0)-LWR, i.e., first select the characteristic parameters of 
the highly enriched elements of the MOX fuel, then determine the number of 
MOX rods of each module, and finally determine the number of fresh refueling 
modules. 


In order to avoid the formation of power peaks on the boundary between the 
MOX elements and the UOz elements, generally there is a buffer zone in the 
low-enrichment MOX region. 


The MOX element has a melting point 20°C lower than that of the U0, element, 
and its thermal expansion is approximately 1 percent larger; in addition, 
its y self-screening effect is better and its heat production is higher. 
These properties partially compensate for the poor thermal properties of the 
MOX fuel. 


The release of gas from fission products and the expansion of U0» or MOX 
particles are not affected y adding Pu to U, but they are affected by the 
microcrystalline structurel!], 


When the power exceeds 430 W/cm, the plutonium in the fuel rod begins to 

shift toward the high temperature region; but this is not a problem in MOX-LWR 
because its linear power level generally is less than 430 W/cm; also, the 
self-screening effect tends to lower the temperature at the center of the 
MOX-LWR. 


Based on its 20 years of experience, the BN Co. in Belgium found that: 1) 
MOX fuel contained in a stainless-steel shell performs better than v0, !41], 

2) there is less mechanical interaction between the MOX fuel plates and the 
shell than UO? fuel elements; 3) MOX has better creep characteristics than 
U0? because of the smaller plate-shell interaction; 4) MOX has slightly 

lower production of 135Xe than U0. 135Xe is an extremely strong neutron 
absorber, but the effect on reactivity is quite small because MOX has larger 
neutron absorption cross section than UO2, thus LWR's using MOX elements 
operate more steadily; 5) MOX elements cause the moderator temperature 
coefficient to become more negative, which may result in a slight loss of the 
reflectivity margin when the reactor is shut down [5], 6) because of the more 
negative Doppler coefficient and moderator temperature coefficient of the MOX 
elements, the power of MOX-LWR decreases more slowly than U07-LWR; this leads 
to better self-screening and better self-modulation; 7) the y scan spectrum 
of discarded elements showed evidence of good power eo when the 
burnout of “plutonium island" elements is very high ; 


A study done by CEA of France showed that in a mixed loading of MOX-U0O, 
elements, the neutron flux is discontinuous, therefore, the reactor physics 
method of calculation should be modified [7], 


The research done by Sweden showed that the slow neutrons emitted by 
plutonium isotopes are two-thirds less than those emitted by uranium isotopes. 
Therefore, MOX elements will respond faster to a change in reactivity than U0) 
elements. This would have an undesirable effect in case of an accident 

such as ejection of control rod in PWR or dropping of control rod in BWR, 

but because of the large negative Doppler coefficient, the accident will 
terminate very quickly [8], 
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Fig. 3 presents a reactivity equilibrium diagram of PWR proposed by the West 
German Power Plant Union (KWU); it shows a comparison of reactivity in a PWR 
when different fuel elements are used [5]. 


Fig. 3 Reactivity Equilibrium Diagram Dec ea 1P Bap. 
for Different Fuel Modules in a PWR 0% SSOnTot 








Key: 
(1) emergency shut-down Doppler effect 
(2) reactivity due to decreasing 
coolant temperature 
(3) returning to critical state 
(4) reactor shut-down margin 
(5) jammed control rod 
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III. Fuel Cycle 


The fuel cycle is generally divided into the front and rear segments. The 
front segment includes mining, mineral »rocessing, enrichment and element 
fabrication; the rear segment includes temporary storage of nuclear residues, 
post-processing, conversion, waste treatment and permanent disposal. The 
recycling of U and Pu in LWR reduces the load on mining, processing and 
enrichment of natural uranium, but it poses some new problems for element 
fabrication, post-processing and conversion. This section briefly discusses 
the difference between MOX and UO? and the resulting effect on fuel cycle 
industrial units. 


1. Post-Processing 


The main differences between MOX and UO? elements with the same burnout and 
the same cooling period are: a) MOX elements have different levels of fission 
products, approximately 30-70 percent higher in 106Ru, !06ph, 85xr and 1291, 
and approximately 10 percent lower in 95Zr-95Nb; b) MOX elements have 
different levels of the actinium series--3 percent less in U, 1-3 times more 
in Pu, Am and Cm; c) there are also differences in physical-chemical 
properties such as thermal property, heat-conductivity, solubility, melting 
point, and density. These differences have the following effects on post- 
processing: 


(1) Segment Processing. There are no major changes in such procedures as 
transportation, fuel removal, storage, disassembly of modules and cutting 
of fuel rods between MOX-LWR elements and U07-LWR elements. 


The dissolving of MOX waste elements is a major problem. If the 

conventional method of mechanically mixing the U0» and Pu0? powder is used, 
there will be a large particles which are difficult to dissolve even in 
boiling, high-concentration HNO, solution. The best way to solve this problem 
is to choose a fabrication process for MOX such that the U and Pu atoms can 

be interchanged in the lattice. Research showed that the following process 

is quite satisfactory: a) grinding the UO and PuO») particles to a highly 
refined powder (less than 0.5) and carefully mixing them; b) applying the 
conversion process to the mixture. 
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(2) Extraction. The Pu content in MOX is 2-3 times higher than that in UO?, 
but since the distribution ratio of Pu(IV) in the 30 percent TBP-OK is smaller, 
appropriate measures must be taken if the level of Pu loss in the extracted 
liquid waste is to be maintained at the same level of U loss in processing 

UO? fuel. If the number of extraction stages is increased or the extraction 
Stage is lengthened, then high-acid (>3M) operation should be used between the 
extraction stage and the washing stage instead of 2MHNO3, so that Pu(IV) can 
be extracted into the organic phase. High-acid operation will reduce the 
cleaning action for Zr/Nb, but will enhance the cleaning of Ru/Rh; hence 

it is compatible with the variation in fission product spectrum (Zr, Nb 

less than 10 percent, Ru, Rh greater than 70 percent). 


The separation of U and Pu can be accomplished by existing reduction 
techniques, but most of these techniques can be improved, or a combination of 
two different techniques can be used 


The same purification process is used for plutonium as in UO?-LWR, i.e., the 
electrolysis reduction process. But because of the 2-3 times increase in 
plutonium content, the design of post-processing plant for U07-LWR should be 
modified to avoid adverse effect on its processing capability. 





In the purification of uranium, it is required that the total « specific 
radiation of the uranium product be less than 15,000 dpm/gU. Because of the 
increased super uranium content in the wastes of MOX-LWR the purification 
process must be strengthened to improve its cleaning ability. 


(3) The Three Wastes of Post-Processing. Sixty to seventy percent of the 
total tritium is collected on the shell of the MOX-LWR elements, which is 

30 percent more than U0?-LWR. The structural components and the head and tail 
blocks of the fuel modules do not differ from those of UO9-LWR. The amount 

and radiation level of the dissolved residues are approximately the same as 
those of UO»-LWR. The concentration of fission products in the high- 

radiation liquid waste is approximately 20 percent higher than that of UO»>-LWR, 
and the amount of heat released is also higher; this will have an effect on the 
heat released is also higher; this will have an effect on the heat dissipation 
and cooling desing of the HAWC container. When glass solidifies, its 

envelope volume must be reduced. Because of the large plutonium content, the 
amount of medium and low radiation liquid waste will increase, and the amount 
of-plutonium loss in the liquid waste is 1 percent less than that of U07-LWR. 
The amount of escaped volatile nuclear elements is also higher than that of 
UO»-LWR, specifically, 3H is 30 percent higher. 85xkr is 40 percent higher, and 
128] is 30 percent higher. There is an obvious drop in solvent effectiveness 
due to increases in a dosage and in 6/y dosage; therefore, the solvent 
processing technique and equipment must be chosen with care. 


(4) Feasibility of Post-Processing. The above analysis yields the 
following conclusions: a) there is no basic difference between MOX-LWR and 
U02-LWR post-processing except for quantitative variations; b) the post- 
processing of MOX-LWR waste elements can be accomplished by making simple 
modifications to the existing methods for processing UO_LWR. 
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2. Conversion 


The process of transforming Pu(NO3), and U0)(NO3)> obtained from post- 
processing or a mixture of the two into Pu0:, UO> or (Pu,U)0» is referred to 
as the conversion process. It is an essential link between post-processing and 
element refabrication. There are two approaches used to produce (Pu,U)09: 
"conversion followed by mixing" and "mixing followed by conversion". 


(1) “Conversion Followed by Mixing". This approach involves first converting 
uranium and plutonium into MO», then mixing them; it has been in use for 
almost 30 years. In the early days, the MO» particles were large in size and 
poor in uniformity; they were difficult to dissolve even in boiling high- 
concentration HNO3 solution. Only by adding HF to HNO. was it possible to 
reduce the insoluble residues to a level of less than | percent, but the 
introduction of HF leads to severe corrosion of the equipment. Subsequently, 
it was discovered that the above difficulties can be overcome by pulverizing 
the MO? particles to a size smaller than 0.5y so they are in a free flowing 
state. 


(2) “Mixing Followed by Conversion". This approach involves first mixing 
Pu(NO3)4 and U02(NO3)2, then converting the mixed sc lution into (Pu,U)0>. 
This method provides a direct and simple way to dissolve MOX. Currently, 
three such techniques are being studied extensively: the sol-gel technique, 
the direct denitration technique, and the joint precipitation technique. 


(i) The sol-gel technique. Based on the requirement of reactor core design, 
the Pu(NO3)4 and the U02(NO3) 9, solutions are mixed according to a prescribed 
proportion and evaporated in vacuum; then, by adding a certain amount of 
reagent to the mixture and passing the mixture through a pulse jet nozzle, 

“gel balls" are formed. Finally, they are washed with NH,OH, dried, inspected, 
baked, and sintered to produce the raw material for element fabrication. 


(ii) The direct denitration technique. In this case, (U, Pu) NO, is converted 
directly into (U,Pu)O9 in either a vulcanization bed or a microwave oven. 
Both techniques share the following common features: 


a) the mixed U. Pu solution can be converted directly; b) the conversion 
process is simple; c) the product requires further pulverization; d) very 
little liquid wastes are produced. But the vulcanization technique has a 
number of advantages over the microwave technique [12, 13), 


(111i) the Joint precipitation technique. This is an advanced technique for 
producing MOX which uses NH,OH as the precipitation system. The Pu(IV) 
oxalate system is ideal for plutonium conversion, but it is not suited for 
joint preciptiation because the oxalate of Pu precipitates faster than that 
of U, thereby causing changes in the MOX composition. The Alkem Co. of West 
Germany developed the AUPuC process which uses NH4,OH to produce a solid 
mixture is then filtered or directly sintered with the slush in a 
vulcanization bed. 


The precipitated MOX powder, like the product from microwave denitration, are 
small in size (0.67) and have large surface ration (9.5 m2/g); they are also 
highly reactive and can readily be made into tablets. The precipitation 
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technique produces about twice as much ammonium nitrate waste water as the 
sol-gel technique. The joint precipitation technique has been used by some 
countries, and each batch produces approximately 1-2 tons of fuel ; 


3. Element Fabrication 


Most of today's large-scale MCX-LWR element fabrication processes are based on 
the “conversion followed by Mixing" method. The procedures are similar to 
those used for U07 element fabrication, and the equipment are essentially the 
same. 


In recent years, the United States had established a "Consolidated Fuel 
Post-Processing Plan’ (CFRP) under which the Oak Hill Laboratory, the 
Hanford Laboratory and the Pacific Northwest Laboratory are jointly studying 
the “gel ball-tablet" technique, and have already produced high-quality 
elements with high density (90-96% TD), good microcrystalline structure and 
high mechanical strength. 


There are other techniques in which the gel balls from the sol-gel technique 
are sintered and inserted into the fuel tubes through vibration to form the 
element rods. Practical experience with this technique is still quite limited; 
the United States plans to use it as the technique for advanced fuel 
fabrication in the future. 


The large amount of plutonium required for the fabrication of MOX elements 
produces some special problems; specifically, they are: a) the problem of 
self-generated heat of Pu; b) the problem of dissolving of Pu; c) the problem 
of hydrogen content; and d) the problem of contamination of plutonium dust 11}, 
These problems can be overcome by choosing the appropriate procedures, using the 
proper equipment and materials, and through intelligent design and management. 





IV. Safety Analysis 


Safety analysis is a complicated issue which covers a wide range of different 
areas. The U.S. Government and the European Common Entity had organized a 
special group of scientists who had devoted a considerable amount of effort 
in carrying out the safety analyses of recycling U and Pu in LWR and 
publishing an extensive report. 


1. Main Conclusions of the European Common Entity 


The European Common Entity made the assumption that by the end of this century, 
the total nuclear power capacity of the Entity will be 150 GWe, and the 
primary power source will be LWR. It is estimated that at that time, 10 tons 
of fissionable plutonium will be recycled each year, producing 400 tons of MOX 
fuel, which is equivalent to about 10 percent of the LWR fuel. These MOX 

will be recycled in 50 of the 150 power plants; they account for 30 percent 

of the fuel used in these reactors, which will be operating in a self- 
sustaining mode, i.e., they produce all the plutonium being recycled. Based 

on this scenario, the group of European scientists reached the following 
conclusions: a) with regard to public radiation level, it is 51 personsievert 
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with recycling, 70 person-sievert without recycling, and the natural 
background is 400,000 person-sievert, thus the level with recycling is less 
than the level without recycling which is less than the natural background; 

b) with regard to individual radiation level for a professional worker, 

there is a 10 percent with recycling, but the professional group radiation 
level is reduced by 10 person-sievert per 10 ton of plutonium; c) in any case, 
the radiation dosage is far below the allowable value specified in the 

safety standards; d) with regard to non-radioactive environmental effects, 
there is no significant difference between the recycling case and the non- 
recycling case in terms of land use requirement, transportation, non-radioactive 
wastes and heat release. 


2. Main Conclusions of the U.S. Government 


After several years of detailed analyses and calculations, the U.S. NRC 

issued a final environmental report which addresses the issues of health 
protection, safety, and environmental effects. In particular, they 

compared five different fuel cycling modes: model: immediate post-processing 
and U recycling, delayed Pu recycling; mode 2: delayed post-processing followed 
by U,Pu recycling; mode 3: immediate U,Pu recycling; mode 4: U recycling, 

but no Pu recycling; mode 5: no U or Pu recycling. 


The basic assumptions for these estimates are as follows: a) time period: 
1975-2000; b) nuclear power capacity: 5X10!1 W by the year 2000; c) total 
power output: 3.5X10!3 KWrh. 


The key results of the report are summarized in Tables 3, 4 and 5. It is seen 
that: a) the effect of U,Pu recycling in the U.S. on public radiation dosage 
and on the cancer death rate is only 1 percent of the natural background; 

b) the differences in the effect on either public radiation dosage or cancer 
death rate between the different modes are far below the fluctuation levels 

of the natural background. 


Table 3. Effect of U, Pu Recycling on Human Radiation Dosage in the 
United States 


Key: 
1. rad dosage, 106 person rem 6. outside residents 
2. mode 7. entire population 
3. natural background 8. foreign countries 
4. U.S. 9. whole world (including U.S.) 
5. professional workers 


Table 4. Effect of U, Pu Recycling on Human Health in the United States 


Key: 
1. effect on health 5. general public in U.S. 
2. mode 6. whole world (including U.S.) 
3. natural background 7. genetic defects 


4. cancer death rate 
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Table 5. Effect of U, Pu Recycling on the U.S. Environment* 


Key: 
(1) other radioactive elements 
(2) residue from ore dressing 
(3) super uranium solid 
(4) high-radiation solid waste 
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Environmental Mode 1.2.3 Mode § Mode ¢ 
Factors LoS aie Se SP ae > eee 
land area, acres 3.4xie 4.0x 10 5.0x 10" 
water, gallons UX awa 
& giseiratEe 2.9x 10" 2.910" 2.9% 10" 
os coal, tons** a ~ 9.9x 10° 9.0x 10 9.0x 10 
$= gas, 105 BTU taxi axa 
= 3 oil, gallons 2.0% 10 2.0x 10" : 19x10" 
Elec GW-Y ext  gexie .exie 
mRe kato 2.5x 10" 2exie 
tT) ne 7, 1.4xiet 
a erie ene 
oT axis OXE gt 
Pu(a) x eee axe 
Pu(f)  Lax108 Tae toe 3.0x 107 
TrPa at si EO 
3H | 65x10" 6.4x 10" xi 
p> : “C | 1.2« 10° 12x 108 a 4.3 10* 
SS OKs 13x10" 1.3 10° 2.6 10° 
& “ "Sr 1.8 10" 1.810" 2.5x107 
"Tc 4.5 10" 5.3 107 — 
a | 1.1« 10" 11x10" i ; 
wit 4x10 3.3 10°  6oxde 
_ iV uae  u_ 
Factory 6 (2)  $.9xi0 6.9% 10" rs}? 
waste (3) 1.5 108 1.3« 10° | - 
(me) (4) 6.5109 6.5x 10 6.5% 10 
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* These effects include: mining, ore dressing, UF; convers.on, U element 
fabrication, MOX element fabrication, reactors, post-processing, 
transportation, waste management, and waste fuel storage. 


** Includes the coal used by fuel recyling plants and by coal-fired power 
plants which provide 2/3 of the electricity. 


V. Analysis of Social and Economic Benefits 


The economic benefits of recycling U. Pu in LWR are quite obvious; even 
more apparent are the social benefits such as conservation of resources 
and safety considerations. 


1. Savings in Natural Uranium and Separation Power 


Because of the different assumptions used by different countries, the 
estimated amount of savings vary from one country to another. Fig. 6 shows 
that approximately 35-40 percent of natural uranium can be saved by 
recycling U, Pu in LWR. For a 1300 MWe PWR, this is equivalent to 80 tons 
of natural uranium each year. 


2. A Good Way to Solve the Pu Storage Problem 


The purpose of recycling Pu in LWR is to replenish the depleted uranium; 
but once it reaches a certain level, power generation can be maintained by 
the self-generated plutonium, and the plutonium regeneration process will 
become self-sustaining. This form of “energy-producing storage” is clearly 
more desirable than the conventional form of storage. The cost savings 
realized by a 1000 MWe LWR by using recycled Pu over PuO? storage are as 
follows: approximately (1.9-4.8) x10 dollars per year for 5-year storage; 
approximately (3.1-6.2)x106 dollars per year for 15-year storage (based on 
1983 U.S. dollar). These data show that it is uneconomical to store Pu for 
20-40 years until FBR becomes operational. 


Table 6. Savings in Natural Uranium Due to Pu Recycling 
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P. recyc. VU recyc. P.+U recyc. 

a) | | 7 

(INFCE) (1980) 20% 20% —_— 

Ger. (BMFT) (1981.10.) 20—30% 15—20% (— 2%) 0% 
(NRC) (1976.9) - 10% 22% 

0.5. (ERDA) (1976.12) - _ 26% 
Japan “(NAIG) (1982) - 17-19% = 
France (CEA—CEN) (1982) ~ 18% - 
Belgium (BN) (1083) - _ ~35% 
(2) (U.R.G) (817) -_ = 20% 
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Key: 
(1) International Nuclear Fuel Cycle Evaluation 


(2) United Post-Processing Company 
3. Moderate Reduction in Fuel Cost 


It is INFCE's opinion that from a narrow economic view point, the economic 
benefits of recycling U, Pu in LWR are quite limited, but if one considers 
the depleting resources, rising uranium prices, and high cost of developing 
FBR's, recycling U and Pu will become very attractive. 


The conclusion of NRC is: The cost saving of recycling U, Pu over the 
single-pass mode of operation is about 8.1 percent. 


4. Becoming Self Reliant in Nuclear Fuel and Protecting the Fuel Supply 


Recycling U, Pu in LWR also helps realizing the goal of becoming self- 
reliant in nuclear fuel and protecting the fuel supply. This is significant 
for countries which are developing their LWR industry but have very limited 
uranium resources. 


VI. A Quick Review of China's Strategy on LWR Fuel Cycle 


According to China's nuclear power development plan, the near-term efforts 
will concentrate on building million KW-class PWR's and on developing FBR's. 
One of the important issues that affects the operation and economic benefits 
of nuclear power plents is to decide on a fuel cycle strategy. In the 
following sections, we shall present an analysis of this issue based on 
China's current conditions. 


l. Three Different Modes of Fuel Cycle 


During the several decades when LWR's transition to FBR's, there are 
three possible modes of processing UO7-LWR nuclear wastes (Fig. 4). Their 
advantages and disadvantages are summarized below: 


Mode 1: a) After several decades of storing post-processed Pu, all the 24!pPu, 
which is 1/7-1/8 of the fissionable plutonium, will become 24lAm and lost; 

b) Because of the strong «, y, n radiation of 241 am, special production lines 
must be built to meet super-plutonium requirements, and costly purification 
procedures will be required for fabricating FBR elements. 


Mode 2: a) One seventh to one eighth of the fissionable plutonium will be 
lost; b) The problems of storage and heat release, corrosion, monitoring, 
damage and retrieval must be solved. 


Mode 3: a) There is basically no loss in fissionable plutonium; b) A 35-40 
percent saving in natural uranium and separation power can be achieved; c) The 
cost of fuel cycle is reduced by 10 percent; d) It contributes toward self- 
reliance in nuclear fuel and protection of fuel supply; e) The energy of Pu 
can be fully utilized, i.e., 3/5 of the time before FBR becomes operational, 
it is used to generate electricity. 
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The above analysis shows that mode 3 provides the best social and economic 
benefits. 


Fig 4. Three Possible Modes of LWR Fuel Cycle 
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2. Current Status of China's Nuclear Fuel Industry 


(1) Nuclear Resources. China's proven resources that have industrial values 
can provide enough fuel for military use and for the operation of 15 million-KW 
PWR's for 30 years [15]. By around 2020, the uranium resources will be mostly 
depleted [16, 17], Also, the average grade of China's uranium ore is quite 

low and the costs of mining and separation power are high. 


(2) Post Processing Industry. China's post-processing industry has 

almost 20 years of operating experience, and has a team of well-trained 
researchers, designers and operators, who had laid a solid foundation for LWR 
post-processing. 


(3) Plutonium Industry. China has accumulated considerable experience in 
the operation of the plutonium industry to be able to develop an industry of 
MOX element fabrication. 








(4) Economic Conditions. China is still an economically underdeveloped 
country with per capita output of less than 400 dollars per year. Therefore, 
it is not possible for China to spend a great deal of money in this century 
to develop fast reactors or to devote its resources to large-scale development 
of PWR nuclear power plants. Under these circumstances, U and Pu 

recycling is attractive because only a moderate of money is needed for the 
design, fabrication and inspection of MOX-LWR elements. 


Based on the above discussions, the author believes that recycling U and Pu 
in LWR is consistent with China's current resource, economic and industrial 
conditions, and is consistent with the guidelines of “becoming self-reliant 
in nuclear fuel, reducing the cost of power generation"; it is also feasible 
from technical and safety considerations and provides stimulant for the 
development of FBR's. Therefore, it is recommended that recycling U and Pu 
be adopted as China's fuel cycle strategy. 
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SINGLE BOARD DIGITAL CONTROL SYSTEM FOR EPITAXIAL REACTOR 
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[Article by Liang Renqiu [2733 0117 7264] and Huang Shengjun [7806 0524 6511] 
of Qinghua University: "Single Board Direct Digital Control System for Epi- 
taxial Reactor"; manuscript submitted in January 1985, revised manuscript 
received in August 1985] 


[Text] Abstract: This paper discusses the structure and characteristics of a 
single board digital control system for a vapor phase epitaxial reactor, the 
mathematical model established, and the two designs for the digital controller. 
The control software and the operation of the system are also described. The 
system has already been successfully placed on a production line. 


l. Introduction 


In manufacturing modern integrated semiconductor circuits, very stringent 
requirements are imposed on the thickness and dopant levels in the epitzaial 
layer on silicon. For this reason, the temperature of the silicon wafer in 
the SVEN the gas flow rate and the processing time must be accurately con- 
trolled. In other countries, small computers are being used for direct digi- 
tal control.* In 1983, China completed the development of direct digital con- 
trol by using microprocessors in the laboratory. However, it is too costly. 

At the end of 1983, Qinghua University and Beijing Electronic Tube Plant joint- 
ly developed a direct digital control system for a Chinese made epitaxial 
reactor, resulting in a single board control system with some special features. 
The system uses a Model TP-801 single board computer for direct digital con- 
trol ard program control of a vapor phase silicon epitaxial furnace. It was 
brought on line in late 1984 and has been working satisfactorily to date. 


2. Structure and Characteristics of the System 


The structure of the system is shown in Figure 1. The single board computer 
sends out the digital temperature signal. Through the 8 bit D/A converter 
(DAC), the optoelectronic coupling (isolation) circuit and the switch S, the 
three phase ac voltage regulating thyristor device can be controlled. The ac 
output of this device is stepped up, rectified and filtered, and then sent to 
the LC oscillator which generates a 400 kHz current to heat the epitaxial 





furnace. The voltage regulator, step-up transformer, rectifier and filter also 
form a local negative feedback system to regulate the silicon wafer temperature. 
The wafer temperature is measured and displayed by a non-contact optoelectronic 


infrared thermometer. 


The system employs seven isolated switching ports. Sampling cycles and proces- 
sing periods are generated by using 4 channels of the Z80-CTC and some software 
program. In addition to the 6 and 4 digit displays on the single board compu- 
ter and thermometer, respectively, the system has a 11 digit display where 4 
digits are used to show the processing time, 3 for HCl flow rate and 2 each 

for PC13 and SiC14 flow rates. The PIO on the computer is used for displaying 
various epitaxial steps and "flow" platterns for the gases, as well as for an 
optical-acoustic alarm. 
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Figure 1. Structure of Single Board Direct Digital Control 
System for Vapor Phase Epitaxial Reactor. 


Key: 
l. Atmosphere 12. Automatic 
2. Electromagnetic valve 13. Voltage loop 
3. Silicon wafer 14. Digital display 
4. Graphite boat 15. Single board computer 
5. Exhaust 16. Optoelectronic isolation 
6. Silicon epitaxial furnace 17. Voltage regulator 
7. Switching parameter output 18. Rectifier and filter 
8. Optoelectronic isolation 19, Oscillator 
9. Flow display 20. Switching parameter input 
10. Photo-acoustic signal 21. Optoelectronic insulation 
ll, Manual 22. Temperature monitor 
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Based on Figure 1 we know that the system is a serial direct digital control 
system consisting of an analog inner loop and a digital outer loop. It is 
also a program controlled system for various valves. The system has the fol- 
lowing special features: 1) In the neighborhood of 1200°C, the temperature 
control error is less than + 4°C. 2) Through man-machine dialogue, major epi- 
taxial parameters can be changed at any time to quickly change over to another 
type of integrated circuit with ease. 3) Silicon wafer temperature, processing 
durations and gas flow rates are digitally displayed. If the temperature ex- 
ceeds its limits or the computer is out of order, an optical-acoustic alarm 
will be activated so that the system can immediately be switched to manual 
operation. 4) A method involving the use of a bias and program switching is 
used to realize the high resolution conversion of temperature signal by a low 
resolution ADC. 5) The program and some data are loaded in an EPROM; 

6) Effective measures such as low pass filtering, isolation, de-coupling, 
shielding and grounding are used to resist interference. Even in the presence 
of intense interference by the 400 kHz high frequency power supply and spatial 
electric field, it can operate reliably on the production line. 7) The equip- 
ment and instruments used, including the DC-8B microcomputer control device, 
temperature probe and high frequency power supply, are made in China. The 
costs are very low. 


3. Establishment of the Mathematica! Model 


In order to find a transfer for the epitaxial furnace including the inner vol- 
tage loop, the LC oscillator and the temperature measuring device, we created 
a rising curve digital measurement system as shown in Figure 2. Near the 
actual operating temperature, the single board computer sends output with sud- 
den increment or decrement. The relevant rising curve is automatically stored 
in internal memory and then printed out. Based on these data, a curve can be 
plotted and the transfer function can be determined. Through plotting and 
computation, the time constant 1] is 45 seconds, the magnification factor Kj 
is 3.44, and the pure time lag 6 is 1 second. The transfer function can be 
written as: 


W »:CS) =Kje~°*/(2,S+ 1) 
=3,44e7*/(45S + 1) (1) 


Figure 2. Digital Measurement System for 
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Because the pure time lag is very small compared to the time constant, the 
system can be treated as a first order inertia system whose transfer function 
is: 








W »,(S) =3.44/(445S +) (2) 


This is the transfer function measured in channel 1 of the ADC. Because the 0 
channel used for PID modulation is also used to amplify the analog voltage to 
be converted by a factor of 4.4 to raise its resolution, the transfer function 
for this temperature control channel should be: 


W (S)=K{ x 4.4/(+ 1) =K,/(4S + 215/458 +1) (3) 


After considering the transfer function Wh(S) of the zero order holder which 
is formed by the DCA input register, the final complete transfer function is: 


W p(S)=W pS) - W,(S)=CK,/(4S+D] - 0Cl—e77*)/S) 
= 15¢1—e77*)/S(45S +1) (4) 
T is the sampling cycle in the above equation which is chosen to be 1 second. 


By performing z transformation for equation (4), the pulse transfer function 
is obtained. 


W »(2) =0.329727*/ (1 —0.9782z7") (5) 


If conditions permit, we may use the least square method to conduct off-line 
recognition to directly dete -mine this pulse transfer function. 


4. Digital Controller 





Based on the requirement in the epitaxial process, the silicon wafer temperature 
to be controlled by the system should follow the pattern shown in Figure 3. From 
room temperature To to point a, the system uses an open loop control to raise 
the temperature at a slope 4; without measuring temperature. From point a to 

b, the temperature is checked in an open loop control to increase the tempera- 
ture at a slope 49. Once the silicon wafer temperature is raised to b (temper- 
ature Tp), the system enters the PID modulating mode. It automatically switches 
channel of the ADC to channel O and begins PDI closed loop control. From point 
c to d, it is the constant temperature segment. It is required to be under PID 
control to ensure the required accuracy in silicon wafer temperature. After the 
epitaxial procedure is completed, i.e. from point d toe and e to f, open loop 
control is used to lower the temperature at slopes 43 and 44, respectively. 





ay Figure 3 
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There are many ways to design the digital temperature controller in segments 
be and cd. Because the system is a definitive single input/single output 
system, it is all right to use the conventional PID method or some practical 
engineering algorithm with pure lag. For instance, it is possible to use the 
following direct digital design method.4»> 


Figure 4 shows the block diagram of the structure of the system. D(z) is the 
pulse transfer function of the digital controller. Y(z), U(z), E(z) and C(z) 
are the z transformation of the system output, given value, control error, and 


digital controller output, respectively. The closed-loop transfer function of 
the system is: 


W .(z)=D(2)W,,(2)/T1+D@)W ,()) (6) 


From equation (6) we can obtain the pulse transfer function of the digital 
controller: 


D(z) =C(2)/E(2) =W..(2)/W,()[1-W (2) (7) 


When the static error of the system is 0, we get the following based on the 
final value theorem: 


e(co) = Lim{(1—2~)U (2) C1 —W.. (2) ]} =0 (8) 
z-1 


Let us assume that the input is a unit step function. The unit step pulse 
function is: 


U(z)=1/Ci-2"") (9) 


To make equation (8) equal to 0, the closed loop pulse transfer function of the 
system should be Wc(z) = z-l. By substituting equations (5) and (9) into (7), 
the specific pulse transfer function of the digital controller is: 


D(z) =C(2)/E(2) =W.. (2) /W, (2) C1 — W. (2) = (3.033 — 2.96627) /Ci—27*) (10) 


This is a PI modulated model. By multiplying equation (10) to C(z)/E(z) and 
then performing an inverse z transformation, we get the following difference 
equation or control scheme which is suited to be computed on a computer. 


e(k) =e(k—1) +3.033e(k) -- 2.966e(k—1) (il) 


This is the theoretical design. We mentioned that the closed loop pulse trans- 
fer function of the system is We(z) = z-l. This indicates that the output of 
the system can vary with respect to the previous input within the same beat 
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(one sampling cycle). Due to the saturation characteristics of the loops in 
the system, this is not attainable. The parameters in equation (11) must be 
corrected accordingly in testing. 
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Figure 4. Block Diagram of the Structure 





As a result of testing, the parameter expressed in equation (11) were found to 
be usable. The temperature variation is less than + 4°C and the transition time 
is under 45 seconds. The parameter was modified to include b2e(k-2) term and 
point differentiation was introduced; i.e. to adopt the following PID modula- 
tion model: 


e(k) =0(k—1) + bge(k) — bye (K— 1) + be CK—2) (12) 


The performance is thus improved. Tests showed that there was no overshoot 
when temperature was raised at a given slope. The static temperature varia- 
tion was less than + 4°C and the transition time was under 35 seconds. When 

a 17°C temperature disturbance was introduced abruptly, the recovery time was 
under 25 seconds. These specifications can satisfy all technical requirements. 


The Dahlin algorithm was also used to design the digital controller®»/, 
Based on this method, when zeroth order register is not taken into account, 
the transfer function is: 


W ,(S) =Kye~**/(, 5+) (13) 


The design is to find an appropriate digital controller D(z) for closed loop 
control. The closed loop transfer function is an inertia loop and a pure lag 
loop in series. In addition, the pure lag time is that of the subject, i.e. 


W CS) =e7* */(tS +1) (14) 


Let the lag time 0 be N times the sampling cycle T, where N is a positive 
integer, i.e. 0 = NT. The pulse transfer function of the closed loop system 
is 
W (2) =20W.(S) + Ci—e“™*)/S) 
=F (fe~*"*/(rS4 1)) - CCl—e7**)/S)} 
=(1—e""*)2-*" 8/1 — 22) (15) 
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The pulse transfer function after taking the zero order register into account 
is: 


W (2) =F (C1 —e77*)/S] + (Kye7***/(4,8 +1) ]} 
= Ky"! (i —e"*/t)/U —e7 2) (18) 


By substituting equations (15) and (16) into equation (7), the pulse transfer 
function of the digital controller can be obtained: 


_C®) _ W ,.() 
D®)=F@) = -W,@li-W.@)) 
o (i—e7**) Cl —e7 7/42) 
~ Kye?) C1 eo 2 — (i eet 


ad —e~*’*)/K, —e~*/**) —( —eT*)eW 2 / KC —e~/**) 
= 1 —e7*/*2-*§~(| —e~*/*)z-*- 











—b,z"* . 
ae (17) 





By cross-multiplying equation (17) with C(z)/E(z), we get 


C(2) (i —a,27* —a,2~-""") = (6,— 6,2") E(®) 


By doing a inverse z-transformation, the algorithm of the digital controller 
can be obtained. 


o(k) =a,e(k—1) +4,0(kK—N —1) + be(K) — bye A—D (18) 


where a,=e7"*, a=l—e~**, by=(1—e7?"")/KiCi—e7"/""), by = bee?” 
(Based on equation (17).) 


When this algorithm was used, some parameters were slightly modified; i.e. 
K1=13, 11=48 sec, T=1 sec, 0=1 sec, N=O/T=1, t=1]/3=16 sec. Based on equa- 
tion (18), the specific control algorithm of the system can be obtained: 


e(k) =0.93040(k —1) +0. 06060(K — 2) +0.22450(K) —0.210900(K—1) (19) 


The result of this algorithm showed that no temperature overshoot at a given 

slope. The static temperature variation was less than + 4°C. The transition 
time was under 45 seconds. The recovery time from an abrupt increase of 50°C 
was under 25 seconds. This method is not sensitive to slight changes in 1], 

Kj and 0. The Dahlin method is a convenient way to design a digital control- 
ler. 
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5. Control Software and Its Operation 


A comprehensive control program, including the control algorithm, switching 
parameters and alarms, was complied for a Z-80 microprocessor. There are 20 
programs including the main program and subroutines. The program takes up 
approximately 4 k in memory and is loaded on an EPROM. 


The operationof the system is shown in Figure 5. The flowchart controlling 
the epitaxial process is shown in Figure 6. The flowchart for the control 
program torealize the temperature pattern shown in Figure 3 is shown in 
Figure 7. 


The system is easy to operate. For example, the power is turned on at the 
beginning of the day. The CPU and PIO are automatically initialized. The 
program and data are loaded in the RAM from the PROM. There are three user 
defined keys as well. By pressing the "set" key, in a man-machine dialogue 
mode, silicon wafer processing temperature, HCl etching time, transition 

layer growth time and epitaxial layer growth time are entered sequentially (in 
decimal numbers). They are displayed again to check for errors. If HCl etch- 
ing or transition layer growth is not required, then just enter 00 minute. By 
pressing the “prepare” key, the lines are "purged." After pressing the "“pro- 
cess" key, epitaxial process begins in furnace number 1. When similar process 
is done later, it is only necessary to press the “process” key besides loading 
and unloading the wafer. The entire epitaxial process will be done automati- 
cally. 
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Figure 5. Operation of the System 
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Figure 6. Flowchart for Epitaxial Process 


4 


1. Press “process” key 9. H2 treatment 

2. Set parameters and clean lines? 10. grow transition layer (or skip it) 
3. no ll. grow epitaxial layer 

4. Return to monitor mode 12. lower temperature 

5. yes 13. cool down 

6. clean reactor 14. issue process finished signal 

7. raise temperature, PID control 15. Return to monitor mode 

8. HCl etching (or skip it) 
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Figure 7. Flowchart of the Temperature Control Program 


Key: 

1. Activate interrupt to protect site 8. greater than TR? 

2. Call sampling subroutine 9. greater than 900°C? 

3. Temperature overshoot? 10. Call PID subroutine 

4. yes ll. Cool down according to y*y-A4 

5. Switch to overshoot treatment 12. raise temp. according to y=4] 

6. rising constant temperature 13. cool down according to y*y-A3 
segment? 14. raise temp. according to y*y+A2 

7. greater than 900°C 15. Shut off interrupt to restore site 


16. Return 
6. Conclusions 


Based on on-site testing and continuous on-line operation, the system performs 
well. It is reliable and meets all the original design requirements. After it 
was put in the production line at Beijing Electronic Tube Plant, the yield of 
epitaxial processes is greatly improved. 
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INTRODUCTION TO PROGRAM VECTORIZATION 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14, No 3, May 86, 
pp 102-109, 113 


[Article by Fan Zhihua [5400 2784 5478] of Changsha Institute of Technology: 
“An Introduction to Vectorization", manuscript received in February 1985] 


[Text] Abstract: A systematic and comprehensive description of the basic con- 
cepts, principles, techniques and applications of the vectorization of serial 
arithmetic is given. 


1. Main Frame Computer and Vector Computer 


Since the 1970's, because of rapid advances in semiconductor and integrated 
circuit technology, a lot of high performance computers emerged. One type of 
computers continues to make breakthroughs in raising the computing speed. The 
average number of arithmetic operations per time unit (second) is increased 
astronomically toward the order of magnitude of billions. This is usually 
called “main frame" computers] (main frame for short). It is not characterized 
by speed alone. 


Main frame computers are the best in terms of performance and size. They are 
the key equipment in solving major problems in scientific computation, engi- 

neering design and data processing. One way to attain super high speed is to 
increase its main frequency. The frequency of the CRAY-1 is 80 MHz.**3 How- 
ever, the primary means is to drastically improve its capability in parallel 

processing. 


The tremendous parallel processing capability of main frame computers is based 
on effective system structures such as pipeline technology, multifunction 
components, sample structure and multiprocessors.**”’ To equip main frames 
with vector mechanisms is the current trend.“ For example, the ILLIAC-IV, 
STAR-100, CYBER-203, CYBER-205, TI-ASC, CRAY-1, CRAY-1S, CRAY X-MP and CRAY~-2S 
developed in the United States, the Model 757 computer? and galaxies computer 
developed in China, and the FACOM VP-100/200 super high speed computer in 
Japan® are vector computers. On the other hand, vector computers are at least 
one to two orders of magnitude higher in parallelism and they are main frame 


computers. 
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2. Vector Mechanism in Vector Computers 


A vector computer is equipped with vector processors. For instance, the FACOM 
VP-100/200 computer develop: ! by Fujitsu in Japan includes vector addition 
operators, vector logic operators, vector multipliers, vector dividers, masked 
registers, load storage and vector registers which are used to protect vector 
operation data (see figure 1).© Based on these hardware resources, a vector 
computer provides its users with vector mechanisms such as vector operation and 
vector chaining. 


A vector is a non-empty ordered set of data of the same type. Vector operation 
deals with vectors. The vector length may be several hundreds, i.e. several 
hundred scalar operations are executed in parallel. From the programming point 
of view, CRAY-1 is equipped with eight vector registers. Each register is made 
of 64 64-bit elements (i.e. vector elements) and has 41 vector commands./~3 
The FACOM VP computer equivalently has 16 times the vector registers as conm- 
pared to the CRAY-l1. It has 82 vector commands .® 


Under the premise of vector operation, vector computers often use the vector 
chaining technique.!»2 When the result of the previous command is the operand 
of the next command (which is frequently so in a program), in the same timing 
cycle, the vector register not only can receive the result from a functional 
component but also can deliver it to another functional component as the source 
for the next operation. Thus, two or more functional components are linked to 
resolve any contradictions related to the operand. 


For example, we are looking for the product of a 5 column matrix X and a 5x125 
matrix Y. The result is a 1x125 matrix Z: 


% ype ™ X eeY ences 


In scalar FORTRAN, two DO loons are used: 


DO 1 / = 1,125 
Z*fy=0 
DO 1 J#1,5 
1 ZU) e ZU) e+ XC YU, D 


Excluding overhead in modifying and decising the loop control variable, it has 
to execute 1254125x5=750 statements. In vector FORTRAN, it only employs a 
simple sequence of six statements: 


Z(1:125) = 0 

Z(12125) = ZC 325) + XC CYC, 12125) 
Z(i: 126) = ZC 125) + XC2) e YC, 12125) 
Z(12126) = Z(1 8125) + XS) © YCS, 15125) 
ZCit 125) = ZC 125) + ACA) © YCA, 15125) 
Z(1: 825) = 2012125) + XG) © VCS, 1125) 
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to be linked for execution. In addition, time overlap is utilized to further 
accelerate the computing speed. 


Parallelism in computer system contains an abundance of contents+’»8. In this 
paper only vector processing is discussed. Based on the classification intro- 
duced by M. J. Flynn in the United States, this parallel mechanism belongs to 
the single instruction multiple data stream system. This hardware based a 
parallel mechanism also requires that the data storage addresses be equally 
spaced. 
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Figure 1. Block Diagram of Processor Part of FACOM VP-100/200 


Key: 

l. Mask register 9. Load etorage 

2. Mask 10. Division 

3. Addition logic operation ll. Maximum 256MB 

4. Vector component 12. Buffer register 
5. Main memory 13. Scalar operator 
6. Load storage 14, Scalar component 
7. Vector storage 15. Channel 

8. Multiplication 16. General register 


Floating point register 
3. Embodiments of Programming Language 


In order to use this hardware based vector mechanism, for the users, the con- 
cern is what descriptive means the programming language can offer.? Since 

the early 1970's, a great deal of studies had been done in vector languages. 
The mainstream is to add vector statements to high level serial languages. 10,11 
For instance, the LRLTRAN on STAR-100, IVTRAN, CFD and GLPYNIR on ILLIAC-IV, 
STARAN on STARAN-IV, BSP FORTRAN on BSP, CRAY-FORTRAN on CRAY~-1, vector FORTRAN 
on Model 757 and Glaxies, and vector ALGOL 60/42+13 are such high level vector 
languages. 
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Because data arrays are the data type used in high level languages, and one- 
dimensional array and its fragments as well as the fragment of a multi- 
dimensional array in a dimension are vectors, array operations are used in the 
expansion. The array processing capability of CRAY FORTRAN makes it possible 
to apply most of the operations defined for scalars in FORTRAN to arrays. 

The four basic arithmetic operations and logic operations can be performed 
among arrays or between an array and a scalar. An array can be assigned. The 
array elements and its operations are conditionally selected. It is used as 
the input parameter of a function or as a process parameter. /° 


In the example in the last section, the statement 
Z(1:125) =0 


is a vector assignment statement. "1:125" is an abbreviation of the three- 
dimensional selection symbol "1:125:1" in which the subscript is selected from 
1 to 125 at an increment of 1. The statement also assigns 0 to Z(1), Z(2), 
and “randomly”. Similarly, the statement: 


Z(12125) = Z(12125) + XC) e YC, 12125) 


is a vector statement which calculates the following 125 multiplications and 
additions and assigns Z(1), Z(2),..., Z(125) in parallel: 


Zi) +XC) YU, YD 
Z(2) +X) * YC, ® 


2125) + XC) © YU, 125) 


The above vector operations are of the single instruction stream type. The 
individual processed data are stored in equally spaced intervals. This is 
ensured by the array compiler and storage distribution rule. 


4. New Theories in Computer Science 


It is known that one mind cannot be focused on two things. The traditional 
human thinking process is serial. Even in parallel events, things are being 
thought one by one. “Taking a walk in the park on Sunday and watching the 
flower show" can be considered as a parallel event. However, when people 
arrange this event, it is still a serial thinking process. The vector computer 
and high level vector language, nevertheless, ask us to describe serial thoughts 
with parallel statements. The following example show that simple substitutions 
will lead to mistakes. 


Let us assume that there is an array A(0:101) in FORTRAN whose initial value 
is A(l)*I1, I*0, 1, ..., 101. Therefore, the serial program 


DO 2/=1, 100 
2 ACI) = AC +1) 





individually assigns A(2) to A(101) to be A(1) to A(100), respectively. The 
corresponding parallel program should be: 


3 AC: 100) = A(2:101) 


The serial program 


DO 4 /=1, 100 
4 AQUI) =ACU- 1D) 


individually assigns A(0) to A(99) as A(1) to A(100). The corresponding 
parallel program should be 


5 A(1: 100) = A(0:98) 


The results of these four programs are shown in Table 1. 


From Table 1 we find out that the vectorization of a pair of similar serial 
programs may end up with totally different results. DO loop 2 and vector 
statement 3 end up with the same results. The results of DO loop 4 and vector 
statement 5 are quite different. 


As another exampie, the double loop 


DO 6 /=1, 100 
DO 6 J=1, 100 
6 BU, J)=BU, J-1)* *2-44# ACU, 1) 


and 
DO 7 J#=i, 100 
DO 7 J#=1, 100 
7 BU, J)=BU, J-1)* ©2-4eACU, J) 


are equivalent. However, the operation of the inner loop of the former program 
requires iterations everywhere. Its parallel algorithm tree is shown in Figure 
2. The depth is 99. The inner loop of the latter program involves independent 
operations. Its parallel algorithm tree is shown in Figure 3. The depth is 1. 
Obviously, the intrinsic parallelism is quite different. 
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Table 1. Results of Four Programs 

























































































Saxe 1.4.1.4 BSeReSPamer 
=F 41 0 | Lanes 4 5 | 6 | 7 | = | 05 9697! 98 | 98 | 100) 101 
Jal i a = = | Ls J 
4 &@ of: ]2]s]a]s5 | 6 | 7 | = | 95] 06) 07) 98 | 90 | 100 | 101 
DO 2 0 2 3 i 4] 5 6 | 7 8 | = | 86 07 | 98 99 | 100 | 101 | 101 
3 ol2e;sial;s}e|}7]8]- 96 | 97 | 98) 99 | 100 | 101 | 101 
DO 4 ofof}o]| e olofot}otl-~ jo 0 | 0} 0 | 0] 0 | 101 
if 
5 0 oj! 2 3 ‘| 5 6 | ~~ (04/95/06) 87 | 98 | 98 | 101 
Key: 
1. Subscript 3. Statement 
2. Contents 4. Initial value 


aaa o computing node 
® number of operations 


Figure 2. Algorithm Tree of Inner Loop of DO 6 Statement 


Figure 3. Algorithm Tree of Inner Loop of DO 7 Statement 


In summary, as vector computers and higher level vector languages become avail- 
able, in order to offer algorithms with higher intrinsic parallelism, the 
development of parallel algorithms is booming. To automatically identify the 
parallelism hidden in serial programs, automatic identification in vector 
operation (vector identification or vectorization) has become a new subject 

in the domain of computer science. 


Parallel algorithms and vector identification go hand in hand. In general, the 
better the parallel algorithm used in the program, the more vector operations 
produced by vector identification will be. An analog can be found in geology 
to describe their relationship. The coal layer was formed in the Permocarboni- 
ferous Period. Today, geological theories are used as weapon in geological 
surveys to search for coal. Parallel algorithm is equivalent to the earth 
environment in the Permocarboniferous Period and vectorization is more or less 
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equivalent to geological survey. The best theory cannot help us locate coal 
in a planet where coal does not exist. Similarly, the best vectorization 
theory cannot find vector operations from programs without any potential. 


5. Basic Vectorization Theory 


If a single instruction single data stream serial computer is considered as a 
dormitory with one seat for a person to dine in, then a vector computer is a 
restaurant with hundreds of seats to serve hundreds of people at the same time. 
The unit is “person * meal", i.e. meal per person. When we receive a reserva- 
tion for 100 person * meals, we are not sure that there are 100 seats simul- 
taneously available in the restaurant. We must first analyze a distribution 
plan. Let us consider the extremes. If there are 1 hundred people eating one 
meal, then we can arrange 100 seats for them to eat it at the same time. If 
one person is to eat 100 meals, then we have to serve one meal at a time even 
though there are empty seats available. It is not even possible to ask him to 
eat two meals together. The reason is obvious. In the former case, 100 meals 
are to be consumed by 100 people. They 2re independent. In the latter case, 
the same person will eat 100 meals. There are rigorous constraints. The meal 
served several hours after breakfast is lunch. 


The vector mechanism in a vector computer his similar characteristics. The 
inner loop of DO 7 in the previous section is equivalent to 100 people eating 
one meal each. The inner loop of DO 6is equivalent to one person eating 100 
meals. The former can be vectorized while the latter cannot be. 


Thus, the purpose of vectorization theory is to automatically identify paral- 
lelism in the source program based on an in depth analysis of the complex 
dependence of the data and to rewrite it in the vector form to obtain a vec- 
torized target program. Obviously, the objective is to improve its efficiency. 
The equivalence of the source and the target program is a necessary premise. 

If vector identification is recognized as the identification work alone, then 
vectorization is more generalized, including the entire task of identification 
and rewriting. 


Based on programming language theory, 16 a source program is statements in the 
source language and a target program is statements in the target language. The 
descriptive capabilities of the source language and target language determine 
the type and complexity of data dependence in the source and target programs; 
which also determine the abiiity of the vectorization theory to identify and 

to rewrite. To date, all vectorization methods are based on FORTRAN, !/~-30 
Most recent results! 0 use FORTRAN77 as the source language!4 and use 
extended FORTRAN 77 with array processing eggenessty such as that of the CRAY 
computer (VFORTRAN) as the target language. + (See Figure 4). The vectoriza- 
tion method thus derived is in general adaptable to various languages. 
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vectorization of program 
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(identify and rewrite) 
Figure 4. Program Transformation in Vectorization 


A review of the history of the study on vectorization shows three areas of in- 
terest: explicit data dependence, implicit data dependence and identification 
instruction. 


As mentioned before, the hardware based vector mechanism has two special fea- 
tures: single instruction multiple data stream and evenly spaced multiple 
data storage. Our programming design experience tells us that the statements 
in a DO loop can easily satisfy the requirements. Each statement in the DO 
loop gives an instruction (or a series of instructions) to be repeatedly 
applied to many data points, forming a single instruction multiple data stream. 
So long as the step remains unchanged, with the aid of the compiling system, 
the data storage happens to be equally spaced. On top of that, the loop con- 
centrates the computing load of the program. Thus, all studies are focused on 
the vectorization of DO loops. Of course, it is also possible to form a loop 
with other statements (such as GOTO). However, based on the principle of pro- 
gramming structure design, 16,31, 32] they can basically be transformed into DO 
loops by the source program itself. 


The data dependence which is caused by the statements in the loop and covers a 
domain limited by the loop is an explicit dependence. Since the 1970's, there 
are many methods to identify explicit data dependence.!/-25, 28-30 They will 
be reviewed in the next section. 


The data dependence which is caused by the statements in the loop and covers a 
domain beyond the loop is an implicit dependence. 28 To make sure that the 
source program is equivalent to the target program, vectorization theory must 
take implicit data dependence into consideration. 


The current vectorization theories do not run the source program in the 
identification stage. Instead, they identify and rewrite in a “static mode.’ 
How to utilize dynamic information in static processing is an issue of great 
concern. Identification guidance instruction is an important way for the user 
to provide auxiliary information, particularly dynamic information which is 
lacking in static identification. 


6. Review of Capability to Identify Explicit Data Dependence 

The subscript tracking method," ->* shed is established on the basis of using 
array elements as analytical nodes to solve the vectorization of all three 
basic modes (link, select and repeat) 16,31, 32 in program structure design. 
Other methods such as coordinate method,!/ hyperplane method!8 and corelation 
analysis method2® use the subscript as the analytical node. 
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To date, it is limited to the domain of the assigned statement loop. Through 
typical examples, their capabilities are reviewed in this section. 


(1) For example, the following is a rigorous loop of assigned statements: 


DO 8i=i, N 
ACI) = BU) * CC) 
CU) =BU-1) 

8 BUI) =ACU+1)°¢D 


All methods described above can be used to identify and rewrite it. 


(2) A loop of assigned statements which can be vectorized along a slant line, 
such as the following!8 


DO 9/=2, 19 
DO9J=2, 2 
9 AU, J)=(AU 41, J+ ACU, J++ AU 1, J+ AG, J-dd 0 02 


can be identified and rewritten by the hyperplane acthod.** 


9 
| 


(3) A loop of assigned statements with a sequence change point, ; such as 


the following 
DO i0 j/=1, N 
Ale 1) = BC) 
10 BBCi) = ACI) +1 


20-22 
a 


can be identified and rewritten by the subscript tracking method nd 


corelation analysis method. ¢ 


(4) A loop of assigned statements whose absolute step is greater than 1 and 
the initial DO variable cannot be estimated,*2»23 such as the following 


M= -Nee2 
DO 11 J/=M, 0, 2 
AQ) =1 
A(2¢1)=2 

i] A(3el)=3 


> 
can be identified and rewritten by the subscript tracking method. ’ 


(5) A loop of assigned statements with dynamic assigned step," such as the 


following 
S = ABS(M - K) 


DO 12 /=1, N, Ss 
A@2¢1) 2 BCI) 
12 Bl) =AU)+1 


can be identified and rewritten by the subscript tracking method. 23 
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(6) A loop of assigued statements with unestimable initial value, final value, 
step and subscript expression coefficient,?? such as the following 


INTEGER M,N, L, S$ 


DO 13 J=M, N, See3 
A(L eI) = BU~L) © C(3e1+1) 
18 BU-L)=CU)e 03-A(Le!) 


can be identified for vectorization in the original form by the subscript 
tracking method. 29 


(7) A loop with logical IF statements or IF blocks,24 such as the following 


DO 14 J=1, 10 
ACI) = BUI) CU) + DU) « «2 
IF (AQ).GT.4) THEN 
A(2¢1+9)=ECI)e ©2+3e ACI) +5 
A(2e1+10)=2¢ FU) © ©2+AU)+1 
G(3 0 1-2) = B20 1-1) © CU) + A201 +9) - AQ eT +10) 
ELSE ;' 
G(3e 1-1) = AU) * AC eI +8) 
G01) = AC) * AQ ei +10) 

i4 END IF 


can be identified and rewritten by the subscript tracking method. 


(8) A loop with a triple control switch, such as the following 


DO 15 /=1, 10 
AC) = BCD) *CU) + DU) « #2 
IFCACI)=4) 18, 17, 16 

16 + A(2*1+9) = EC) © ©2438 ACI) +5 
GOTO 20 

17 = A(201+10)=2e FCI) © © 2+ AC) +1 
GOTO 10 

18 = G(3e 1-2) = B(2 © f= 1) © CUI) + A(2 0 1 +9) = Ao oT +10) 

19 G(3*1=1)= AC) © AQ eT +9) 

20 G(3e1) = ACI) © A(2 01 +10) 

16 CONTINUE 


can be identified and rewritten by the subscript tracking method. *" 
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(9) A loop with unestimable initial value, final value, step and subscript 
coefficient and with logic IF statements and IF blocks,2° such as the following 


INTEGER M,N, L, S 


DO 21 /=-2°M, N+L, S* 2 
IF ((BU-L)+ A(LeI1)).EQ.C(2*I-1)) THEN 
ACL #1) =BU=-L)* ACL *1)+DU +1) 
ELSE 
BUI -L)=CU1) « ©3-A(Le]) 
END IF 
21 DU +1)=B8Ci-L) 


can be identified and quasi-vectorized by the subscript tracking method. >” 


7. Applications of Vectorization Theory 


To realize vectorization theory in engineering applications, the vector identi- 
fier is a key component to improve the operating efficiency of the vector 
computer. A vector identifier is a new software system of considerable size. 
On the vector computer developed by China, the earliest vector identifier is 
still in operation. It is being used in the following areas: 


a. Transform Existing Serial Programs 


Before vector computers became available, people already wrote a large number 
of serial programs which were successfully tested after investing an enormous 
amount of manpower and money. Now, we obviously would like to use these pro- 
grams still “without making any changes." If we do that, since we did not in- 
clude the vector mechanism when we wrote and tested these programs, the 
efficiency will in general be very low, leading to significant waste. Under 
the circumstance, vector identifier plays the role as a piece of “transplant 
software" to alter these serial programs so that they are automatically modi- 
fied into new programs suitable for the operating environment of vector 
computers. 


b. Modify Future Serial Programs 


There are still a variety of computers in operation in parallel to vector compu~- 
ters. Because these computers can only execute serial programs, in order to 
fully utilize the results obtained on these computers, vector computers must 
have the capability to accept new serial programs and run these programs very 
efficiently. This capability comes primarily from the vector identifier. In 
addition, there will be considerable vector computer users who write and test 
their programs on other computers. The programs are run on vector computers 
only after they are compiled and debugged. These users will also introduce a 
lot of new serial programs. The vector identifier also plays the role of soft- 
ware transformation for these programs. 
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Figure 5. Program Vectorization by VFORTRAN 
c. Modify Future Hybrid Programs 


After vector computers are available, many users will write programs and run 
them directly on vector computers. They are not limited by serial FORTRAN. 
The system allows them to use direct vector statements. 


However, as we mentioned before, our traditional way of thinking is serial. 
Unless in extremely simple cases, the most natural and dependable way for a 
programmer to express a parallel operation on a vector computer is still by 
using a serial program. 


Thus, we will be facing a lot of iybrid programs with serial and parallel 
statements. In simple and obvious cases, vector operation will be used. In 
other areas, scalar operations are still being used. For this portion of the 
program, the vector identifier is again an optimization software. It can 
optimize the scalar operations which are masked by the complicated data 
dependencies into vector operations. 


Because FORTRAN/7 is a subset of VFORTRAN, the source program for the former 
must be that of the latter as well. From the viewpoint of program transforma- 
tion, the function of the vector identifier is to extend the vectorization of 
FORTRAN77 to VFORTRAN (Figure 4) to the transformation of VFORTRAN itself 
(Figure 5). 


d. Partially Modify Program Structure 

As mentioned above, vector recognization is somewhat equivalent to a geological 
survey which basically identifies existing parallel components, instead of 
creating parallel components. 

However, there are exceptions. For example, a vector identifier can alter the 
program structure within a small range to vectorize the program in Section 4 
by changing DO loop 6 to 0 loop 7. 


8. Conclusions 


As vector computers become more popular, vectorization studies have both theo- 
retical and practical significance. The essence of vectorization theory Lies 
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in the static recognition of dynamic data dependencies. Its progress not only 
requires knowledge in mathematics and computer science but also relies on the 
vector identifier. It is expected to be widely used in computer science and 
mathematics such as in parallel algorithm, cyclic optimization and data stream 
analysis. 


Vectorization is being perfected and will continue to be an important part of 
computer science. 
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CLONING OF THE P7 5 PROMOTER FROM VACCINIA VIRUS 


Shanghai SHENGWUHUAXUE YU SHENGWUWULI XUIBAO [ACTA BIOCHIMICA ET BIOPHYSICA 
SINICA] in Chinese Vol 18, No 1, Jan 86 pp 74-79 


[Article by Feng Zongming [7458 1350 6900], Wu Xiangfu [0702 4382 3940], Wang 
Yuan [3076 0997], Chu Meijin [0328 5019 3866], and Li Zaiping [2621 6528 1627], 
Shanghai Institute of Biochemistry, Chinese Academy of Science] 


[Text] Abstract: This article reports on the gene segment and cloning of 7.5K 
early protein P7_ 5 promoter from vaccinia virus isolated from the genomic DNA 
of vaccinia virus, strain WR, collected in China. Use of restrictive endo- 
nuclease mapping, Southern hybridization, and DNA sequence analysis has con- 
firmed that the P7 5 promoter has been obtained. Comparison with data on the 
WR strain as reported in foreign literature shows its DNA sequence varies by 

a difference in one nucleotide. As one of the stronger promoters found in 
vaccinia virus, the P7_5 promoter can be very useful for activating foreign 
gene expression. 


The host environment for vaccinia virus used in genetic engineering studies is 
extensive and safe, and the capacity of the vaccinia virus genome to carry for- 
eign genes is also great, thereby permitting insertion of many foreign genes. 
Hopefully, the vaccinia virus can be used this way to produce multivalent vac- 
cine. For this reason, research during the past 2 years on the use of vaccinia 
virus to study the expression of foreign genes has been on the increase. 1-6 

In this approach, expression of the inserted foreign genes need to be control- 
led by the vaccinia virus genome's own promoter gene, and not by any foreign 
promoter. 


The basic purpose of developing such studies is to obtain a strong promoter 
from the vaccinia virus and to allow the expression of foreign genes which is 
controlled by this promoter in the vaccinia virological system to produce 
needed products. This article reports primarily on the cloning and assay of a 
comparatively strong P7,_5 promoter fragment from vaccinia virus, WR strain, col- 
lected in China. It also reports on its nucleotide sequence and its difference 
from known WR strain P7 5 promoter. 
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Materials and Methodology 


1. Vaccinia virus, WR strain. Provided by colleague Yan Zilin [09)7 1311 2651] 
of the Bureau of Pharmaceuticals and Biologicals Assay and Certification, 
Ministry of Health. 


2. Restrictive Endonucleases and T4DNA ligase. Supplied by the Enzyme Tooling 
Section of our laboratory and from Biolabs. 


3. Alpha->“P-daTP ( 3000 Ci/mmol). From Amershan. 


4. DNA isolation. Isolation of vaccinia virus DNA based on method used by 
Esposito et al.’ Isolation and purification of DNA particles, their in vitro 
recombination and transformation, and selection of recombinants etc., are based 
on previous work done in this laboratory.® 


5. DNA probe labeling. Reference made to Rigby's method. ? 

6. DNA gel transfer and hybridization. Reference made to Southern's method. 19 
7. DNA sequence analysis. Based on Maxam-Gilbert's method.!! 

Results 

1. Cloning 


We have completed the cloning cf HindIII fragment of the genome of vaccinia 
virus, WR strain (yet to be published report). However, cloning was difficult, 
because the HindIII-A fragment was too large. Moreover, the linear double 
linear strands at the terminal end of the vaccinia virus’ DNA molecule were 
joined to each other by a covalent bond to form a closed ring. Therefore, the 
HindIII-B and HindIII-C fragments at both ends of the DNA molecule were not 
able to be cloned. Coincidentally, the P7_5 promoter we needed was located on 
the HindIII-C fragment of the left terminal of the vaccinia virus genome. If 
the HindIII-C fragment were to be cleaved with Sall enzyme, the P7_5 promoter 
would be found in the SalI-C, fragment./¢ Toward this end, the recombinant 
plasmid containing the Sall-C) fragment must first be obtained, after which the 
P7.5 promoter could be isolated from it. 


In this experiment, we started by using a BglI-C fragment to replace the 
HindIII-C fragment, because the map positions of these two fragments are quite 
similar. Also, HindIII enzyme cleavage of vaccinia virus, WR strain, yields 
15 segments, whereas cleavage with BglI yields only four segments A, B, C, and 
D,23 the distance between C and B, D segments being rather larger and more 
easily isolated. We first isolated the BglI-C segment. Then the Byll-C seg- 
ment and pBR322 were each cleaved with Sall, after which they were placed to- 
gether for recombination and transformation. LB plates containing ampicillin 
and tetracycline were used to screen-select the Ap'Tc® recombinants. DNA 
extracted from these recombinant plasmids then underwent enzyme cleavage 
analysis to select out the required cloning strain (Fig. 1). 
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Figure 1. Diagram showing the cloning of DNA 
fragment containing the P7 5 promoter 


2. Analysis and Selection of Restrictive Enzyme Cleavage 


Results of enzyme cleavage analysis of the WR strain's genome shows four Sall 
cleavage points on the BglI-C fragment, which yields five segments when cut 
with Sall. At the end of the left segment is the single-stranded ring formed 
by covalent bonding of the double-stranded DNA. At one end of the right seg- 
ment is the Sall cleavage point; at the other, the Bgll cleavage point. The 
DNA in these two segments cannot recombine with pBR322 cleaved by Sall. Using 
an antibiotic-treated LB plate to screen and select Ap'Tc® colonies yields 
only SalI-A, Sal-Cj and Sall-Cj cloned strains!* (Fig. 1), of which the A seg- 
ment is 9 kb in length, while C; and C> segments measure only about 0.9 kb. 
Only analysis of recombinants in the 0.9-kb segments is needed to select out 
the required Sall-C; segment. Colonies are cultured in small quantities for 
quick extraction of the DNA. The recombinant plasmids of 0.9-kb length are 
first selected out. After this, enzyme cleavage analysis, selection, and 
determination are made of the inserted fragment, from which the cloned strain 
containing the P7 5 promoter, that is, the Sall-C; segment, is selected out. 


The results of work by Venkatesan at al]? show that the Sall-C> segment con- 
tains three Hpall cleavage points, which are distributed fairly evenly and 
cut the Co segment of 0.9-kb length into very small fragments. However, the 
Sall-C; segment only has one cleavage point which is displaced toward one end 
of the 0.9 kb long C; segment (Fig. 2A), which may be used for selection pur- 
poses to distinguish between the two. As shown in Figure 2B, after the 0.9 kb 
long recombinant DNA has been enzyme cleaved by Sall, and the two recombinant 
plasmids are found to definitely contain insertion segments about 0.9 kb long 
(Fig. 2B, a and c), adding Hpall will show a segment of about 0.8 kb in the 
DNA of one recombinant (Fig. 23, b), while the other strand is cut into 
numerous small fragments (Fig. 2 B, d). Quite possibly, recombinant a in 
Figure 2B may be the cloned strain containing the Sall-C; segment now called 
pVC,, pVCp being the strain containing the Cp segment. After the DNA of pVC; 
has been digested by enzymes Sphl, Pvul, and Hincll, the resulting segments 
(Fig. 2B, e~g) and sites are similar to those reported by Puckett et ai, 14 

in both the number of cleavage points and sites. 
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Figure 2. Enzyme cleavage electropherogram of pVC, 
a. pVC,-Sall; b. pVC)-Sall-Hpall 
c. PVC9-Sall; d. pVC9-Sali-Hpall 
e~g. Addition to pvC; of SphlI, Pvul, 
and Hincll. 


As for the Rsal cleavage point, according to the restrictive enzyme cleavage 
electropherogram as described by Venkatesan et al,/2 the C} segment only 
indicates one Rsal cleavage point, which is also one of the points used by 
Smith et al* to isolate the 275-bp long P7.5 (P2) promoter. We have found 
three Rsal cleavage points on the Sall-C; segment. Figure 3 shows how the 
Sall-C; fragment (a), separated from pVC; is being digsted by HpaJi(b) and 
Rsal(c); that there is a Hpall cleavage point on the Sall-Co fragment that 
shears out a small fragment of about 100 bp (Fig. 3,b); that shearing by Rsal 
results in four fragments, of which one measures 280 bp. What is interesting 
is the synergistic effect of using enzymes Rsal and Hpall together that 
eliminates the approximately 280-bp fragment in the digestion process. These 
results explain that this fragment is very likely to be the Sall-Cy fragment 


that contains the P. 5 promoter. 
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Figure 3. Enzyme cleavage electropherogram of 
Sall-C, fragment showing 


a. Sali-C, fragment;  b. Sall-C, + Hpall; 
c. Sall-C, + Rsal; d. Sall-C, +Rsal + Hpall 





3. DNA Transfer and Hybridization 


After we had obtained the DNA fragment that very likely contained the P7_5 
promoter, we obtained from Dr. Moss a vector plasmid pCGS20 that contained the 
promoter P7 5 from the 7.5K early protein gene of vaccinia virus. Then 32P was 
used to label the approximately 280-bp fragment isolated from Sall-C) in a 
probe to be hybridized with pGS20. The results are noted in Figure 4 where 

the P7 5 promoter in pGS20 was partially sheared out by EcoRI (Fig. 4A, a), 
confirming the fact that the 32p-labeled probe can hybridize with it, while 

the other parts of pGS20 do not respond to hybridization (Fig. 4A-B, a and a‘). 
Because a small fragment of polymerized connectors is located behind the P7_5 
promoter in pGS20, it shows up as being larger than the approximately 280-bp 
fragment we had obtained (Fig. 4, c and c'). This def‘nitely proves the fact 
that we had found the DNA fragment containing the Pas promoter. 


Figure 4. Hybridization of P7_.5 promoter by 
Southern's method. 


A. a. pGS20 + EcoRI; b. pVC; + Sall; 
c. pVC; + Sall + Rsal. 

B. Respective autoradiographs (a, b, c) of 
a, b, and c in A. 
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4. DNA Sequence Analysis 


Mapping analysis of the restrictive endonucleases described above, and the 
results of DNA gel transfer and hybridization according to Southern's method, 
have confirmed the fact that we had obtained the P7_5 promoter. Requirements 
for the promoter we need must have the distinguishing sequence of RNA poly- 
merase, its synthesis sites and transcription initiation codon, and not contcin 
the translation initiation codon ATG from the original code of the protein gene. 
This way, it will help the foreign protein gene inserted behind it to directly 
utilize its own translation initiation codon ATG to begin translation and avoid 
formation of a mixed peptide chain. For this reason, when the promoter is be- 
ing isolated, the site for accurately excising the promoter must be accurately 
placed between the transcription initiation codon and the translation initiation 
codon ATG of its adjacent protein gene. 


Toward this end, we conducted partial DNA sequence analysis of the P7_5 pro- 
moter. Beginning near the translation initiation codon ATG of the 7.5K protein 
and going toward the transcription initiation codon, we analyzed a total of 142 
nucleotides as shown in Figure 5. When the results are compared with those 
published by Venkatesan et al,)> the only difference was found in the 48th 
nucleotide located upstream along the translation initiation codon ATG5" of 

the 7.5K protein. In their results, this site is occupied by C; in the frag- 
ment we obtained, this is occupied by T. No other differences were noted. 

This experiment to isolate the P7_5 promoter has given us reliable data which 
indicates we need only shear off this DNA fragment and insert it into the 
appropriate site on the vaccinia virus, and the foreign gene spliced downstream 
on its 3' side will be expressed./® 


a. AGTTTTTCAC OCATAAATAA TAAATACAAT Av .TAATTITO 
b. AGTTTTTOAC OOATAAATAA TAAATAOAAT AATTAATTTO 


TOGTAAAAGT AGAAAATATA TTOTAATTITA TIGCACGGTA 
TOGTAAAAGT AGAAAATATA TIOCTAATITA TIGOATGGTA 


AGGAAGTAGA ATCATAAAGA ACAGTACTOA ATCAATAGCA 
AGGAAGTAGA ATCATAAAGA ACAGTAOTOA ATCAATAGOA 





ATCATGAAAO AATATATOOT 00 
ATOATGAAAO AATATATOGT OO 


Figure 5. Partial DNA sequence analysis of P. 5 Promoter. 


a. Venkatesan's results; b. Authors’ results. 


Discussion 


Due to the unique characteristics of the vaccinia virus, the use of vaccinia 
virus DNA as vector units to construct recombinant viruses for studying the 
expression of foreign genes has received wide attention. At present, expres- 
sion of such foreign genes as herpesvirus thymidine kinase, ! hepatitis B sur- 
face antigen,» influenza virus hemagglutinin, 3 herpes simplex virus glyco- 
protein,4 chloroamphenicol acetyltransferase (CAT) ,? and plasmodium sporozoite 
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antigen® in vaccinia virus has become realities. To study these modes of 
expression, a vector plasmid which contains a fragment of vaccinia virus DNA 
into which the foreign gene can be inserted must first be constructed. To 
achieve a more pronounced expression of the foreign gene, a stronger promoter 
from a vaccinia virus fragment is inserted ahead of the foreign gene. Among 
several known vaccinia virus promoters, in vitro experiments have found pro- 
moter P7_5 in the 7.5K early protein gene to be stronger./4 Smith et al2 
have used it to successfully express hepatitis B surface antigen. 





Consequently, we have isolated and purified viral DNA from vaccinia virus, WR 
strain, collected in China, conducted endonuclease mapping analysis, and 
established a gene bank. Following this, we have also cloned a DNA fragment 
containing the /.5K early protein gene, and have analyzed and assayed its pro- 
moter. Using enzyme cleavage mapping, Southern hybridization, and partial DNA 
sequence analysis, we have determined the site of the P7_5 promoter. When, 
the fragment is cut di the Rsal point, an approximately 280-bp fragment is 
obtained between the Rsal and Sall points. As the Rsal cleavage point is 
located 18 bp on the 5° side of the 7.5K protein gene's translation initiation 
codon, it does not contain the 7.5K protein's ATG codon, but contains rather, 
the RNA’s polymerase synthesis locus and its transcription initiation codon. 
We can very conveniently excise this promoter, and insert it into the expres- 
sion plasmid's vaccinia virus tk gene, and construct a transcribed expression 
vector that is controlled by the P75 promoter. We have used it to success- 
fully study expression of the hepatitis B virus surface antigen gene, 16 As 
for our DNA sequence analysis studies where the P7_5 promoter discovered by us 
show a T and C difference with already published results by others, the effect 
of this difference on the strength of the promoter can pose an interesting 
problem. 


Experiments using in vitro systems show the P7_5 promoter to be more effective 
than other promoters. However, the strength of these promoters within the 
genome has not been determined. For this reason, after the location and struc- 
ture of the promoters concerned have been clarified, foreign genes can be con- 
structed behind different promoters, undergo in vivo recombination and be 
placed in vaccinia virus genome for comparing the strength of the various pro- 
moters. Such an approach can be quite significant. 


Using the vaccinia virus promoter to study the regulation and control mechanism 
of gene expression is quite meaningful. Recently, Weir et al!’ took a 366-bp 
fragment of the Py) 7 promoter from a 28K late protein gene and inserted it into 
an early protein--thymidine kinase (tk) gene, to study its effect on the 
expression of the CAT gene inserted behind it, and found CAT activity surfacing 
late, appearing only after DNA replication had begun. This explains the fact 
that this promoter only begins activity at this time. If a DNA inhibitor were 
present, this P;,7 would remain inactive, completely unaffected by the early 
gene environment. On such a short DNA fragment are the on and off transcription 
signals, RNA polymerase recognition signals and synthesis loci and the trans- 
cription initiation site. Not only does selection of a simpler and quicker 
method to screen and select promoters have practical application value, it is 
also significant, through comparative analysis, in understanding the regulation 
and control pattern of true gene expression. 
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Applied Mathematics 
FINITE ELEMENT ANALYSIS FOR TEMPERATURE FIELD PRODUCED BY MOVING HEAT SOURCE 


Chongqing YINGYONG SHUXUE HE LIXUE [APPLIED MATHEMATICS AND MECHANICS] in 
Chinese Vol 7 No 5, May 86 pp 383-400 


[English abstract of article by S. N. Atluri of Georgia Institute of Tech- 
nology; and Kuang Zhenbang [0562 7201 6721] of Xi'an Jiaotong University] 


[Text] Using the moving mesh finite element method, the authors discuss the 
temperature field produced by a moving heat source with variable thermal 
conductivity and with radioactive and convective boundary conditions in a 
wide velocity range. The temperature-time relationships at various velocities 
in the static and moving coordinate systems are studied. The steady-state 
temperature distributions at various velocities in the moving coordinate 
systems are given. The temperature field produced by the plastic deformation 
at the process region (a region very near the crack tip) is also studied, 

and the results show that the highest temperature at the process region is 
lower than 1000°C or 1832°F. (Paper received 9 May 1985.) 
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MORE GENERALIZED HYBRID VARIATIONAL PRINCIPLE AND CORRESPONDING FINITE 
ELEMENT MODEL 


Choneaing YINGYONG SHUXUE HE LIXUE [APPLIED MATHEMATICS AND MECHANICS] in 
Chinese Vol 7 No 5, May 86 pp 443-449 


[English abstract of article by Chen Wanji [7115 8001 0679] of Dalian 
Institute of Technology] 


[Text] According to recent studies of the generalized variational principle 
by Professor Qian Weizhang [6929 0251 7022], a more generalized hybrid varia- 
tional principle for the finite element method is given. Thus, a new kind of 
generalized hybrid element model is established. Using the thin plate bending 
element with varying thickness as an example, various hybrid elements based 


on different generalized variational principles are compared. (Paper 
received 9 April 1985.) 
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Electronics 


NEW LADDER-TYPE ACTIVE BANDPASS FILTER BY USE OF FDNR EMBEDDING, NIC 


Beijing DIANZI KEXUE XUEKAN [JOURNAL OF ELECTRONICS] in Chinese Vol 7, No 6, 
Nov 85 pp 401-410 


[English abstract of article by Li Wenzhe [2621 2429 0772], and Hu Yun [5170 4596] 
of Beijing Institute of Posts and Telecommunications] 


[Text] A new method of realizing elliptic ladder-type active bandpass filters 
is proposed. Comparing with the methods of the same category, it has the 
advantage of using fewer operational amplifiers, so the power consumption can 
be reduced. It also has the advantage of ease in adjustment. The low sensi- 
tivity property of ladder networks is maintained. The stability and the error 
of the fundamental circuits are analysed in detail. A sixth-order bandpass 
circuit is realized with this method and it shows that the experimental results 
agree pretty weli with the theoretical design. (Paper received 3 May 84, 
finalized 10 Aug 84.) 
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THEORETICAL ANALYSIS OF NEW TYPE OF RF SENSOR ON SATELLITE 


Beijing DIANZI KEXUE XUEKAN [JOURNAL OF ELECTRONICS] in Chinese Vol 7, No 6, 
Nov 85 pp 411-420 


[English abstract of article by Zhang Jinbiao [1728 6855 2871] of Tianjin 
Institute of Technology] 


[Text] A theoretical analysis of the performances of RF sensor proposed by 
B.K. Watson et al. (1981) on the large Satellite is carried out in this paper. 
By use of the method given here for dealing with the error signals, it makes 
the RF sensing system be indeed a big performance system. Except the case of 
b=1, y=0 or 1/2, the RF sensor is insensitive to depolarization of the coming 
waves, and because the horizontal channel does not couple with the vertical 
one, so it is obvious that the precision of dynamic tracking of the RF sensor 
on satellite is high and in ground station this RF sensing system will be a 
good one for tracking high-speed vehicles. (Paper received 25 Mar 85, 
finalized 16 Jul 85.) 
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IMPEDANCE MATCHING ELEMENT OF DIELECTRIC IMAGE LINE 


Beijing DIANZI KEXUE XUEKAN [JOURNAL OF ELECTRONICS] in Chinese Vol 7, No 6, 
Nov 85 pp 421-427 


{English abstract of article by Shi Meiqi [0670 5019 3825] of Department of 
Radio Electronics, Beijing University] 


[Text] In this paper the wave impedance of dielectric image line as a function 
of geometry dimensions and frequency is analyzed. Both the necessity and the 
possibility of impedance matching are studied. An optimal matching element is 
put forward. Its bandwidth with VSWRE1.10 can be as wide as 10.5 GHz, when it 
is connected between two dielectric image lines with an impedance ratio of 4./. 
Its properties are obviously better than those of classical type. (Paper 
received 16 May 84, finalized 21 May 85.) 
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NEW FORMULAS FOR DESIGNING FREQUENCY DEVIATION OF CRYSTAL FM OSCILLATOR 


Beijing DIANZI KEXUE XUEKAN [JOURNAL OF ELECTRONICS] in Chinese Vol 7, No 6, 
Nov 85 pp 443-449 


[English abstract of article by Ji Jian [0679 0256] of Nanjing Jinling Zhiye 
University] 


[Text] On the basis of analyzing the oscillatory frequency of crystal oscil- 

lator, the general expression of the frequency deviation of crystal frequency- 
modulation oscillator is derived. The design formulas of the frequency devia- 
tion under optimum load condition and under condition of the modulating stage 

with series negative current feedback are also given. (Paper received 

23 Apr 84, finalized 20 May 85.) 
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I-V CHARACTERISTIC ANALYSIS OF SCHOTTKY AMORPHOUS SILICON SOLAR CELLS 


Beijing DIANZI KEXUE XUEKAN [JOURNAL OF ELECTRONICS] in Chinese Vol 7, No 6, 
Nov 85 pp 458-465 


[English abstract of article by Xu Le [1776 2867] of Department of Physics, 
Nankai University] 


[Text] The present paper introduces an experimental method for measuring the 
width of illuminated and short-circuited Ni/a-Si:H Schottky barriers. The 
current-voltage curves for the Schottky barrier solar cells under AMl, 

100 uW/cm* illumination are calculated by using the parameters determined by 
experiments. The diffusion length of holes in a-Si:H obtained from the 
illuminated I-V curve is consistent with the results measured by the author 
with the surface photovoltage method in 1983. The factors affecting the fill 
factor are analysed on the basis of the calculated results. A comparison of 
the calculated results to the experiment's reveals that the very low fill fac- 
tor of the solar cells measured is due to series and shunt resistances rather 
than the low diffusion length of the holes. (Paper received 16 Apr 84, 
finalized 20 May 85.) 
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CIRCULARLY POLARIZED BACKFIRE ANTENNA WITH HIGH GAIN 


Beijing DIANZI KEXUE XUEKAN [JOURNAL OF ELECTRONICS] in Chinese Vol 7, No 6, 
Nov 85 pp 473-476 


{English abstract of article by Song Ximing [1345 6932 2494) of The Huanghe 
Machine-Building Factory] 


(Text) In this paper, a new type of high-gain circular polarized backfire 
antenna is presented. It is developed on the basis of equivalent reflector 
method proposed by the author (1979, 1981). A circular polarized wave is ob- 
tained by use of two orthogonal dipoles. The lengths of the two dipoles are 
so chosen that the real parts of their input admittances differ in 90°. The 
transverse dimension of the antenna is 2.24A, the axial dimension is 2A, the 
weight of the complete antenna assembly is only 1.2 kg. The gain measured is 
17.5 dB, the axial ratio is 0.5 dB, VSWR is 1.18 at centre frequency. (Paper 
received 23 Apr 84, finalized 8 May 85.) 
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Metallurgy 


FORMATION OF (y+y') EUTECTIC AND CONTROL OF o-PHASE IN HIGH Al-Ti CAST Ni-BASE 
SUPERALLOY 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
18 Apr 86 pp A93-Al100 


[English abstract of article by Zhu Yaoxiao [2612 5069 1366], et al., of the 
Institute of Metal Research, Chinese Academy of Sciences, Shenyang] 


[Text] An investigation has been made of the various phases formed in order 
aud composed under different temperatures during the process of soliditication 
and segregation of alloying elements in a high Al-Ti cast Ni-base superalloy. 
The formation of (y+y') eutectic and the control of o-phase pre :ipitation are 
discussed. During solidification of the alloy, the dendritic segregation of 
Al and V is almost unespied, and that of Ti, Cr and Mo is positive, with that 
vf Ti being the most serious. This may impel the formation of a segregation 
zone at the eutectic reaction and the enrichment with the o-phase forming 
elements, e.g., Cr, Mo and Co, in the remaining liquid. Therefore, the 
increasing of the o-phase transformation may be tended. However, the o-phase 
may be avoided either by adjusting the ratio Al/Ti, taking advantage of the 
difference of the segregation behavior of Al and Ti, or by controlling the 
cooling rate after casting. (Paper received 21 November 1983; finalized 

15 May 1985.) 
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RELATIONSHIP BETWEEN y+c TRANSFORMATION TEMPERATURE M., AND COMPOSITION OF 
METASTABLE AUSTENITE REGION IN Fe-Mn-Al-Cr SYSTEM 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
18 Apr 86 pp Al101-Al108 


{English abstract of article by Tian Xing [3944 5281], et al., of the Institute 
of Metal Research, Chinese Academy of Sciences, Shenyang] 


[Text] The y*e martensite transformation temperature Ms of 40 metastable 
austenitic alloys in an Fe-Mn-Al-Cr system was determined by dilatometry. 

As a result of the multicomponent linear regression, a formula for a linear 
relationship between Meg and the alloy composition has been developed in that 
each wt-percent of Al, C and Cr decreases 51, 339 and 8°C of M.s of the Fe-Mn 
alloy respectively. This is in agreement with the results of previous allied 
works and Ishida's calculation of the influence of alloying elements on the 
driving force of the y*e transformation in Fe-Mn alloys. A brief discussion 
on the effects of Al, C and Cr on M.s is also made which refers to the 
mechanism of stacking fault for the y*e martensite transformation. (Paper 
received 18 October 1984; finalized 4 January 1985.) 


REFERENCES 


1 Eichelman, G.H.; Hull, F.C. Trans, ASM, 45 (1953), 77. 

2 Monkman, F.C.; Cuff, F.B.; Grant, N.J. Met. Progr., 7 (1957), 94. 

3 Pickering, F.B. Physical Metallurgy and the Design of Steels, Appl. Sci. Publ. London, 
1978, p. 228. 

4 £E4, PHAR 1979— 19 SLL TAT. BR, 7.221, HRS LPS, SST Hew 
tt, 1984. 

S$ HBS, GAA, SKFR. 7 (1964), 2S. 

6 SRERFTR. KAFRATH. 

7 &£S3%, FBR, SRR, 19 (1983), A253. 

8 £B+, SUFR, 19 (1983), A262. 

9 Troiano, A.R.; Mcguire, F.T., Trans. ASM, 31 (1943), 340. 

10 Holden, A.; Bolton, J.D; Petty, E.R. J.Jron Steel inst., 209 (1971), 721. 

ll Bx. RAAFRARH. 

12 Handbook on Materials for Superconducting Machinery, Metal and Ceramic Information Center, 
Publ. 1974, p. 5. 1. 1. 

13 Cina, B. Acta Metall, 6 (1958), 748. 

14 Ishida, K.; Nishizawa, T, Trans. Jpn Inst. Met., 18 (1974), 225. 

iS Andrews, K. W., J.Jron Steel Inst., 203 (1965), 721. 

16 Kato, T; Fukui, S.; Fujikura, M.; Ishida, K, Trans. Iron Steel Inst. Jpn, 16 (1976), 673. 

17 Yoshimura, H.; Yamada, N.; Yada, H.; Honma, H.; Ito, T, Troms. Irom Steel Inst. Jpn, 16 
(1976). 98. 

18 White, C.H.; Honeycombe, R.W.K., J. /ron Steel Inst, 200 (1962), 457. 

19 Gras £47: ER, Sum, L®, MPwRt. 1964, Pp. 241. 

20 Breedis, J.F. Trans. Metall. Soc, AIME, 230 (1964), 1583. 

21 Breedis, J. F.; Kaufman, L., Metall, Trans, 2 (1971), 2359. 

22 Ishida, K. Scr. Metall,, 11 (1977), 237. 


94 




















EFFECT OF HYDROGEN ON CREEP IN AUSTENITIC STAINLESS STEEL 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
18 Apr 86 pp A115-Al118 


[English abstract of article by Chu Wuyang [5969 2976 2254], et al., of 
Beijing University of Iron and Steel Technology] 


[Text] Dynamic charging does not influence the creep rate of type 321 
austenitic stainless steel if the initial strain is less than 2 percent 

in which the amounts of evolving hydrogen after the test were below 1.5 ppm 
(by wt). However, the creep rate increased markedly when the initial strain 
was larger than 8 percent or when the specimen was precharged for a long 

time. In this case, the corresponding amount of evolving hydrogen was 5 to 12 
ppm (by wt). (Paper received 1 February 1985.) 
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BEHAVIOR OF Hf IN SOLIDIFICATION OF CAST Ni-BASE SUPERALLOYS 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 


18 Apr 86 pp Al119-A1l24 


[English abstract of article by Zheng Yunrong [6774 6663 2837] of the Institute 


of Aeronautical Materials, Beijing] 


[Text } 


tion. 


alloys. 
form of Hf-rich phases. 


eutectic y'. 


Results of quantitative metallography, SEM, EDAX and EMPA studies of 
Hf on the solidification of cast Ni-base superalloys are presented. 
is found to lower the liquid temperature of the superalloys, with the final 
solidification temperature even dropping to 1130°C. 
the range between the temperature lost during the interdendritic capillary 
feeding action and that of the solid, as well as decreasing the liquid 
content necessary for linking the interdendritic pools in the late solidifica- 
The Hf-rich melts have superior fluidity, wettability and skin effects. 
These are the factors in increasing the castability of Hf-containing super- 
The Hf is not found in the dendritic arms of the superalloys, but is 
exceedingly concentrated in the narrow interdendritic zones, usually in the 
The Hf content in primary Hf-rich phases may be 
arranged in decreasing order of MC(2), (Hf, Ti), SC, Ni<Hf, MC(,) and 

The Hf may depress the formation temperature of primary phases. 
The solidification temperatures of all Hf-rich phases are lower than 1250°C. 
These phases may precipitate from the molten pools with Hf concentration over 
13 wt-percent at late solidification and may have an intergrowth feature. 


(Paper received 21 September 1984; finalized 27 March 1985.) 
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ON RAPIDLY SOLIDIFIED Cu-Zn-Al SHAPE MEMORY ALLOY 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
18 Apr 86 pp Al25-Al129 


[English abstract of article by Jin Jialing [6855 0857 7117], et al., of 
Shanghai Institute of Iron and Steel Research] 


[Text] The shape memory properties, microstructure and influence of heat 
treatment on them of Cu-24.9 wt-percent Zn-4.5 wt-percent Al alloy ribbon 
of different sizes prepared by the rapid solidification technique have been 
investigated by means of o-c and R-T curves, optical metallography, X-ray 
diffraction, SEM and TEM. The results show that the alloy ribbon is of a 
fine-grained 8 structure, about 20 ym originally, and grows to about 100 ym 
after water quenching from 720°C for 10 min. As for the alloy with a §+M 
structure at room temperature, its optimum shape memory properties may be 
obtained and the embrittlement by grain coarsening may be avoided. The 
complete memory recovery may be realized with strain of no more than 8 percent. 
The formation of granular relief and martensite on the specimen surface seems 
to be reversible during deformation under loading. (Paper received 

6 October 1984; finalized 12 March 1985.) 
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EQUILIBRIA OF Ce-Al-O AND Nd-Al-O IN MOLTEN LRON 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA) in Chinese Vol 22 No 2, 
18 Apr 86 pp A139-Al148 


{English abstract of article by Diao Shusheng [0431 3219 3932], et al., of 
Reijing University of Iron and Steel Technology] 


[Text] The equilibrium of the Ce-Al-O or Nd-Al-O systems in molten iron was 
studied by separately smelting the radioactive isotope of ‘*'Ce or '*’Nd in 
an Al,O, crucible filled with pure iron to form ***CeAl0, or **’NdAlO, around 
its inner wall. The content of dissolved Ce or Nd in iron was determined by 
means of radioassay and electrolysis in an organic electrolyte. The analysis 
of the dissolved Al in iron was obtained by colorimetry, followed by calcula- 
tions. The activity of oxygen in a dissolved state in liquid iron was 
directly measured by solid electrolyte sensors made of Zr0,(MgO) tube. y 
extrapolation of the data obtained, the temperature dependence of the 
equilibrium constant of RE-Al-O in molten iron, Kppajo;, or of the intevaction 
coefficient, eft, may be described as: 


For the reaction CeAl0;,(.) = [Ce] + [Al] + 3[0] 
lg Kcealos = ~37900/T + 5.84 
AGRoa1Q 04= 173400 - 26.7T cal/mol 


eAl = 20500/T - 13.53 
For reaction NdAlU,(,) = [Nd] + [Al] + 3[0) 


lg Kndalo; * ~24700/T + 0.69 
AG dAlO, = 113030 - 3.16T cal/mol 


Al . - 
Na 13500/T 9.33 


(Paper received 14 July 1984; finalized 20 January 1985.) 
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PRECIPITATION DIAGRAM OF [Ce]-[0]-({S] IN LIQUID IRON AND ITS APPLICATION 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
18 Apr 86 pp Al149-Al155 


[English abstract of article by Dong Yuanchi [5516 0337 4654], et al., of 
Beijing University of Iron and Steel Technology] 


[Text] A space precipitation diagram of [Ce]-[O0]-[S] in pure liquid iron has 
been constructed by the proper selection of recent and correct experimental 
thermodynamic data collected from allied equilibrium systems of RE-O, RE-S 

and RE-O-S. The quantitative correlation between the activity of oxygen or 
sulphur and the constitution of RE inclusions formed has been evaluated 
theoretically by thermodynamic calculations. This was verified well by the 
experimental results from which the correlation between the initial activities 
of oxygen and sulphur in liquid iron and the constitution of RE inclusions 
formed was also determined. (Paper received 17 April 1984; finalized 

13 January 1985.) 
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EQUILIBRIUM OF Fe-La-O-S SYSTEM 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2 
18 Apr 86 pp A156-Al62 


[English abstract of article by Chen Dong [7115 0392], et al., of Be i jing 
University of Iron and Steel Technology] 


[Text] A study was made of the equilibrium of the Fe-La-O-S system. The 
oxygen activity was directly measured by the solid electrolyte sensors made of 
a Zr0,(MgO) tube. The stability of La,0)S in the electrolytic separation 
and the content of La as a dissolved state in liquid iron were examined by 
means of radioassay sad electrolytic isolation of inclusions in organic 
electrolytes under low temperatures. By extrapolation of the data obtained 
the equilibrium constant of La,0,S may be described as: 


J 


1550°C KLa202S = 3.81 x 107+? eka = -4,] 
1600°C Kia.02s * 6.45 x 107}? ek@ = -4.8 
1650°C Kia.0,s = 7-41 x 107'9 ek@ = -4.2 


La(,) = [La] 
AGS). = -136700 + 59.82T, Yj, = 0.33 (1873 K) 


It seems that the deoxysulphurizing effect of La is more intense than that of 
Ce, Nd and Y. (Paper received 18 June 1984; finalized 15 December 1984.) 
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ESTIMATE OF FRACTURE TOUGHNESS FOR POLYCOMPONENT COMPOSITE MATERIALS 


Beijing JINSHU XUEBAO | ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
i8 Apr 86 pp Al75-Al180 


[English abstract of article by Li Zhonghua 621 0022 5478] of Xi'an 
Jiaotong University] 


[Text] The expressions for estimating fracture 


toughness of the polycomponent 
composite materials have been derived 


im terms ot the fracture toughnesses and 
volume fractions of each individual component included. The expression, 
conforming to the mixed rule, of the upper (m-~) or lower (n=1) bound value 

is also derived based on limit analysis. (Paper received 18 Jun- 1984; 
finalized 8 November 1984. ) 
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KINETICS OF COMBUSTION OF PULVERIZED COAL IN RACEWAY BEFORE TUYERE OF BLAST 
FURNACE 


Beijing JINSHU XUEBAO (ACTA METALLURCICA SINICA] in Chinese Vol 22 No 2, 
18 Ap: 86 pp B49-B62 


[English abstract of article by Yang Yongyi [2799 3057 1355], et al., of 
Beijing University of Iron and Steel Technology! 


{Text! A combustion efficiency study was made of partially burned pulverized 
coal, sampled from the raceway before the tuyere of two blast furnaces. Re- 
sults show that the injected coal started to volatilize and burn just as it 
entered the blowpipe, which seems to be an important preheating and precom- 
bustion space before the raceway; however, it could not be burned completely 
in front of the tuyere. In order to replace more coke by injecting more coal, 
the laboratory simulation test on the combustion kinetics of the pulverized 
coal used by the two blast furnaces under conditions prior to tuyeres indicated 
it was technically possible by: pulverizing the coal down to 40 ym; enriching 
the blast with oxygen up to 30-40 percent; further raising the blast tempera~ 
ture; improving the injected coal distribution early in the blowpipe; con- 
Sistently injecting the coal at all tuyeres and selecting a variety of coal 
which would be more combustible, e.g., bituminous coal is superior to anthra- 
cite. Petrographic observations of coal samples prior to injection and after 
being partially burned revealed that the bituminous coal burned violently and 
Simultaneously, not only on the surcace but also in the interior of its 
particles, while the anthracite burned mostly on the surface. Based on the 
theory of heat-transfer and mass-transfer, a set of mathematical models for 
the combustion of pulverized coal has been derived. By substituting the 
combustion conditions of both blast furnaces and the laboratory test furnace 
into the model and with the aid of a computer, a close approximation to 

what happens in actual practice has been obtained. (Paper received 

25 February 1985; finalized 27 June 1985.) 
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LIQUIDUS AND INTERMETALLIC COMPOUND IN Al-RICH RECION OF Al-Mg-Ce SYSTEM 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
18 Apr 86 pp B63-B67 


[English abstract of article by Zheng Chaogui [6774 2600 6311], et al., of 
the Department of Chemistry, Beijing University] 


[Text] The liquidus of the Al-rich corner in the Al-Mg-Ce phase diagram and 
the intermetallic compound, if any, formed in this region have been studied 
by thermal, EMPA, chemical analysis, metallography and X-ray examination. 

An intermetallic compound in the soiid phase was observed to form a solid 
solution of a rather extensive range. The pseudo-ternary system Al-Mg<Al,- 
CeAl, is of the ternary eutectic type. Three surfaces corresponding to the 
primary crystallization of a-Al, MgsAl, and CeAl, respectively, three 
umivariant curves associated with secondary crystallization and a ternary 
eutectic with a composition of about 12.6 wt-percent Ce and 26.0 wt-percent Mg 
at 446°C were detected. (Paper received 20 August 1984; finalized 

l February 1985.) 
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COMPUTERIZED SIMULATION OF STRUCTURE OF CRYOLITE MELT 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
18 Apr 86 pp B68-B74 


\English abstract of article by Xu Chi [1776 7459], et al., of Shanghai 
Institute of Metallurgy, Chinese Academy of Sciences; Li Qingzhi [2621 

1987 5347], et al., of Shanghai Institute of Computational Technology; 

Shen Shiying [3088 2514 5391] of Northeast Institute of Technology, Shenyang] 


(Text] The structure of molten cryolite has been investigated by computerized 
Simulation using the Monte Carlo method. The radial distribution functions 
and some parameters describing the local structure of 3NaF-ALF, melt at 1283 K 
have been calculated. It is concluded that a part of the F~ ions in ALFj~ 
moves away from Al** ions after the melting of cryolite, forming xNa*-yF~ 
clusters with the Na* ions. The removal of a part of the F~ from ALFi~ makes 
it disintegrate into some local structures, such as AIF; or AIFi~. Some Al’* 
ions are connected to each other by Al-F-Al bridges, forming ionic clusters, 
e.g. AloF> , etc. It is also concluded that the distribution of the free 
space in molten cryolite is not uniform. There are more fissures and holes 

in xNa*-yF~ clusters and between AIF, and xNa*-yF~ clusters. (Paper received 
7? November 1984.) 
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ON THE CRYSTALLIZATION OF AMORPHOUS Ni-BASE ALLOY 


Beijing JINSHU KUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
18 Apr 86 pp 375-879 


[English abstract of article by Wang Yuming [3769 3558 2494], et al., of 
Jilin University, Changchun} 


[Text] The crystallization of amorphous Ni,g,Cr;Fe;Si,B,; alloy was investigated 
under dynamic heating and isothermal annealing conditions. Identical crystal- 
lization and unaltered phases were revealed under both treatments. However, 
the initial crystallization was mainly due to the nucleation of crystals uader 
dynamic heating, while crystal growth played a major role under isothermal 
annealing. It was also found that the microhardness approached its highest 
value just before the specimen was totally crystallized. (Paper received 

5 July 1984; finalized 19 January 1985.) 
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CORRELATION BETWEEN STRUCTURE OF ELECTROCHEMICALLY DEPOSITED (Fe, Co, Ni)-P 
LAYERS AND PROCESS CONDITIONS 


Beijing JINSHU XUEBAO [ACTA METALLURGICA SINICA] in Chinese Vol 22 No 2, 
18 Apr 86 pp B80-B89 


[English abstract of article by Jiang Xiaoxia [1203 2556 7209!, et al., of 
the Institute of Metal Research, Chinese Academy of Sciences, Shenyang] 


[Text] A study has been made of the structure »f electrochemically deposited 
amorphous (Fe, Co, Ni)-P layers in relation to the process conditions, e.g., 
pH value, current density, metallic salt concentration, etc., by means of 
X-ray analysis, SEM and EPMA. Results show that the effect of the process on 
the structure is realized by changing the P content of the layer. When the 

P contained is over 8 percent, below 3 percent or between 3 and 8 peccent, a 
layer may be constructed as an amorphous alloy, supersaturated solution or 
mixed structure respectively. The major factors affecting the structure of 
the deposited layer are found to be the concentration of NaH,PO, and pH value 
of the bath. The structure of the layer may also be influenced by the poten- 
tial of electro-deposition. The Fe-P, Co-P or Ni-P system is similar in 
amorphous structure when the highest thickness of the atomic face stack is 
2.0, 2.5 or 6.5 nm respectively. (Paper received 6 January 1984; finalized 
15 October 1984.) 
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Nuclear Physics 


NEJTRON EMISSION BEHAVIOR IN LIGHT-CHARGED-PARTICLE ACCOMPANIED FISSION OF 
252c¢ 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 289-296 


[English abstract of article by Han Hongyin [7281 3163 6892], et al., of the 
Institute of Atomic Energy, Beijing] 


[Text] The reutron emission probabilities for long range alpha (LRA), 
helium-3, triton and proton accompanied fission of *°*Cf are measured in a 
three-parameter experiment in which a liquid scintillation detector and a 
telescope are used to record the number of prompt neutrons and the energy of 
light charged particles. The average number of neutrons per fission for 
various particle accompanied fission events mentioned above are 3.13 + 0.02, 
3.09 + 0.09, 2.95 + 0.05 and 3.24 + 0.07 respectively. The differentials of 
the average number of neutrons with respect to the energy of LRA particles 
and tritons are ~0.037 + 0.003 and -0.039 + 0.008 neutron/MeV, respectively. 
The dependence of the width of neutron emission probability on LRA energy, 
described by Diven's empirical formula, and the correlations between the 
neutron number and triton energy have not yet been reported in the 
literature. 
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°&,°0,°2Ni(q,p) THREE-NUCLEON TRANSFER REACTIONS AND q OPTICAL POTENTIAL 
AMBIGUITIES 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 297-306 


[English abstract of article by Wang Yuanda [3769 6678 1129] and Bao Xiumin 
[7637 4423 2404], et al., of the Institute of Atomic Energy, Beijing] 


[Text] The differential cross sections are measured using a 26.0 MeV a particle 
for °°.®*Nif(a,a) °%.°*Ni and °*,*®*Ni(a,p) ®'>*®°Cu reactions as well as a 

25.4 MeV a particle for ®°Ni(a,a)®°Ni and ®°Ni(a,p)®*Cu reactions. Consistent 
calculations with the optical model and ZR DWBA are made for (a,a) and (a,p) 
reactions by using single, two, three and four nucleon optical potential 
parameters. For elastic scattering due to the a optical potential ambiguities, 
all the above optical potentials can reproduce the experimental angular 
distributions. However, the single, two and three nucleon potential, including 
Baird's mass systematics and Chang's energy systematics of a potentials, 
obviously cannot provide a reasonable fitting with the (a,p) reaction experi- 
mental data. Only the results from the four nucleon potential are in good 
agreement with the (a,p) reaction experimental data. This reveals that in 

the a-particle induced transfer reactions, the real depth of the a-nucleus 
optical potential should be rather deep. 
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CERIUM CROSS SECTION BETWEEN 12 AND 18 MeV 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 307-312 


[English abstract of article by Teng Dan [3326 0030], et al., of the Institute 
of Atomic Energy, Beijing] 


[Text] The cross sections for ***Ce(n,2n)***Ce and **°Ce(n,2n)!?%9Ce from 

12 to 18 MeV, '?*Ce(n,2n)*?2°Ce and '*°Ce(n,p)'*°La from 13 to 15 MeV, and 
138Ce(n,2n)*?7’™Ce at 14.07 MeV are measured by the activation method. The 
neutron fluence is determined by the associated particle counting at 

14.59 MeV for the '*?,!*°Ce(n,1In) reactions. The '?*,'*®Ce(n,2n) and '*°Ce(n,p) 
reactions are determined relative to the cross sections of the '**»!*®Ce(n,2n) 
reactions at 14.07 MeV. The measured values are 1981, 1766, 925, 1244 and 

6.70 mb for the '*?2,/%°,238,138Ce(n 2n) and **°Ce(n,p) reactions respectively. 
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MICROSCOPIC ANALYSIS OF (p,n) AND (p,p) SCATTERING ON TIN ISOTOPES 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 313-319 


[English abstract of article by Ning Pingzhi [1337 1627 3112] of Nankai 
University, Tianjin] 


[Text] Microscopic calculations are performed for the (p,n) charge exchange 
reaction and (p,p) elastic scattering on '!*Sn and '?"Sn at a proton energy 
of 24.5 MeV. The optical potentials and the isospin-coupling potential are 
generated from a density~dependent complex effective nucleon-nucleon inter- 
action and the nucleon density distributions of the target nucleus. The 
(p,n) and (p,p) differential cross sections are calculated simultaneously 

by solving the coupled Lane equations exactly. The results of these micro- 
scopic calculations yield reasonable agreement with the measured values, even 
though they contain no adjustable parameters. 
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VALIDITY OF M-3Y FORCE EQUIVALENT G-MATRIX ELEMENT FOR CALCULATIONS OF 
NUCLEAR STRUCTURE IN s-d SHEI.L 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 320-324 


[English abstract of article by Song Honggiu [1345 1347 4428], et al., of 
Shanghai Institute of Nuclear Research, Chinese Academy of Sciences] 


[Text] The matrix elements of the M-3Y force are adopted as the equivalent 
G-matrix elements and the Folded Diagram Method is used to calcuiate the 
spectra of '°0 and *®F. The results show that the matrix elements of the 
M-3Y force as the equivalent G-matrix elements are suitable for microscopic 
calculations of the nuclei in the s-d shell. 
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TRANSITION POTENTIAL OF QUARK-ANTIQUARK PAIR (MESON) EXCHANGE IN CONSTITUENT 
QUARK MODEL 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 325-331 


[English abstract of article by He Hanxin [0149 3352 2450] of the Institute 
of Atomic Energy, Beijing] 


[Text] The amplitude to the lowest order for the exchange of quark-antiquark 
pairs (q-q) between two three-quark clusters (nucleons) are obtained by QCD 
perturbation theory when the two nucleons overlap, and thus the transition 
potential for the exchange of the color-singlet (q-G) pairs (mesons) between 
the two nucleons is defined. The non-relativistic approximation of this 
transition potential is also discussed and the expression of the potential 

is consistent with the one-boson exchange one. The nucleon-meson vertex 
structure is consistently included in this potential due to the quark struc- 
ture of the hadron. 
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SHELL MODEL CALCULATIONS OF °°Zr WITH WOODS-SAXON POTENTIAL 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 332-337 


[English abstract of article by Yuan Haiji [5913 3189 7535] of Hangzhou 
University; and Shen Qingbiao [3947 1987 1753] of The Institute of Atomic 
Energy, Beijing] 


[Text] In this paper a set of local Woods-Saxon potential parameters for the 
asymmetric nucleus *°Zr are adjusted to perform the single-particle shell 
model calculations. The calculated results of the single-particle energies 
and the root mean square radius of the charge distribution are in good agree- 
ment with the experimental data. Based on previous studies, the nucleon 
densities, kinetic energy densities and spin densities of the asymmetric 
nucleus °°Zr are calculated. Their shape and values are presented in this 
paper. 
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DEFORMED HF STATES OF EVEN TITANIUM AND CHROMIUM 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 338-343 


[English abstract of article by Liao Jizhi [1675 4949 1807] of Sichuan 
University] 


[Text] Using KB matrix elements, the deformed Hartree-Fock calculations for 
#4, %6,%8,5°7T] and *®,°°Cr nuclei are performed. The oblate, prolate and 
triaxial deformed solutions, as well as the solutions of some particle-hole 
excited configurations, are obtained. The results show that there is a form 
transition from mass number 48 to 50, and that the single-particle energy 
spectra are different not only for different mass numbers of the nuclei, but 
also for different configurations. 
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GLOBAL ANALYSIS OF NEUTRON DIFFERENTIAL ELASTIC CROSS SECTION CALCULATIONS 
WITH MICROSCOPIC OPTICAL POTENTIAL 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 344-351 


[English abstract of article by Tian Ye [3944 6851], et al., of the Institute 
of Atomic Energy, Beijing; and Liu Wei [0491 3555], et al., of Northwest 
University, Xi‘ an] 


[Text] Neutron differential elastic cross sections are calculated by the 
microscopic optical potential (S-MOP) and compared with available and indepen- 
dent experimental data as well as those calculated with Greenlee's phenomeno- 
logical potential. The global analysis shows that for the original Skyrme 
forces, such as GS2 and SKa, without adjusting any parameters, S-MOP can 
reproduce the tendency of angular distributions, positions and absolute values 
of peaks and valleys of experimental data in the neutron energy region of 

1-26 MeV for nuclei ranging from '*C to **°Pu. Thus, it is shown that S-MOP 
has a capability for predicting neutron differential elastic cross sections in 
optical model calculations for a nucleon even-even nucleus in a rather wide 
range of energies and nuclei. 
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FIFTH-ORDER ANGULAR ABERRATIONS OF INHOMOGENEOUS MAGNETIC FIELD WITH AXIAL 
SYMMETRY AND THEIR CORRECTION 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 352-360 


[English abstract of article by Lu Hongyou [0712 3163 3731], et al., of the 
Institute of Atomic Energy, Beijing! 


[Text] The fifth-order approximate solution is derived from the equation of 
motion for a charged particle in inhomogeneous magnetic field with axial 
symmetry by the method of successive approximation. The angular aberrations 
of different orders and the coefficients of the magnetic field to correct 
these angular aberrations are given. 
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TRANSPORT OF BEAM PHASE SPACES WITH NONCENTROSYMMETRIC CONFIGURATION THROUGH 
LINEAR FIELD 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 361-368 


[English abstract of article by Fu Shinian [0265 0013 1628], et al., of the 
institute of Atomic Energy, Beijing] 


[Text] The transport theory of the beam phase spaces with noncentrosymmetric 
polygonal configuration x + bjx' - aj = 0 (i = 1, 2, ..., mn) as well as non- 
centrosymmetric smooth convex close curve configuration f(x, x') = 0 through 
the linear field is developed. The formulas of the beam envelope of such 
beam phase spaces through the linear field are derived. Detailed calculations 
are worked out for typical examples to show the practical application of the 
theory. 
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SCHEME FOR RAISING BEAM BLOW-UP (BBU) THRESHOLD CURRENT OF CONSTANT 
IMPEDANCE STRUCTURE 


Beijing YUANZIHE WULI [CHINESE JOURNAL OF NUCLEAR PHYSICS] in Chinese Vol 7 
No 4, Nov 85 pp 369-378 


[English abstract of article by Yao Chongguo [1202 0339 0948] of Nanjing 
University] 


[Text] This paper presents a design scheme for a constant impedance structure 
that raises the BBU threshold current but does not degrade the performance of 
the structure. Two holes opened on disks in certain positions on the structure 
can stop the EH,, wave that brings about BBU, but they have little effect on 
the TM,, wave that accelerates electrons. Their function is equivalent to 
dividing the whole structure into r segments (n is usually two to four, 
depending on the length of the structure), so that the BBU threshold current 
roughly increases by the same factor of n. 


9717 
CSO: 4009/4 


118 














Nuclear Radiochemistry 


SYNERGISTIC EXTRACTION OF U, Th, Nd, Y¥Y BY BINARY SYSTEM OF AB SPECIES 
COMPOSED OF PMBP-TOPO--A STUDY USING TITRATION CALORIMETRY, KARL FISCHER 
TITRATION AND SOLVENT EXTRACTION 


Beijing HE HUAXUE YU FANCSHE HUAXUE [JOURNAL OF NUCLEAR AND RADIOCHEMISTRY] 
in Chinese Vol 8 No 1, Feb 86 pp 1-7 


[English abstract of article by Wang Wenging [3769 2429 3237], et al., of the 
Department of Technical Physics, Beijing University; Li Xingfu [2621 6821 
1133], et al., of the Institute of High Energy Physics, Chinese Academy of 
Sciences, Beijing; and Sun Pengnian [1327 7720 1628] of the University of 
Science and Technology of China, Hefei] 


[Text] The present study is a detailed investigation of the synergistic 
extraction of uranium, thorium, neodymium and yttrium with a binary system 

of 1-phenyl-3-methy1-4-benzoyl-pyrazolone-5(PMBP) and trioctyl-phosphine oxide 
(TOPO) using a direct calorimetric titration method together with Karl Fischer 
titration and solvent extraction. 


Calorimetric titration and temperature dependence studies of the extraction 
method allow calculation of the enthalpy and entropy changes associated with 
the synergic reaction. From the magnitude and signs of the AS values, it is 
possible to determine whether synergism is associated with the replacement 
of solvate water by the adduct molecule TOPO with no change in the coordina- 
tion number, or with an expansion of the coordination sphere upon addition 
of the adduct molecule but no replacement of water. Uranium, thorium, 
neodymium and yttrium have been shown to form U0,(PMBP),(TOPO), Th(PMBP),(TOPO), 
Nd(PMBP),(TOPO) and Y(PMBP),(TOPO) respectively. Water analysis is made by 
Karl Fischer titration. The equilibrium constants, enthalpy, entropy and 
free energy changes for the addition reaction are determined. The data 
indicate that the main driving force of the synergistic effect is the 
substitution of the coordinated water molecule by TOPO. 
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STUDY OF PREPARATION OF URANIUM (111) BY ELECTROLYSIS AND KINETICS OF 
REACTION OF U** WITH H* 


Beijing HE HUAXUE YU FANCSHE HUAXUE [JOURNAL OF NUCLEAR AND RADIOCHEMISTRY] 
in Chinese Vol 8 No 1, Feb 86 pp 26-31 


[English abstract of article by Wu Jizong [0702 4949 1250] and Luo Wenzong 
[5012 2429 1350] of the Institute of Atomic Energy, Beijing] 


[Text] Im studying the factors affecting the percentage and rate of reduction 
of tetravalent uranium to the trivale.t state by controlled potential 
electrolysis, it is found that the main factor is the cell constant 8, which 
depends on the design of the cell system. The bigger the 8, the higher the 
percentage and rate of the reduction. The kinetics of the redox reaction of 
u** ions with H* ions in hydrochloric acid determined by controlled potential 
electrolysis follows a second order reaction and can be represented as: 


-d(u?+)/de = k{u**)[H*]. 
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INVESTIGATION OF STABILITY OF U30, 


Beijing HE HUAXUE YU FANGSHE HUAXUE [JOURNAL OF NUCLEAR AND RADIOCHEMISTRY]} 
in Chinese Vol 8 No 1, Feb 86 pp 53-55, 7 


[English abstract of article by Feng Guoning [1409 0948 1337], et al., of 
the Uranium Mining Research Institute, Hengyang] 


[Text] The source of the weight increase of U;0, during storage is 
investigated. The main source comes from the absorption of water in the 
atmosphere. Probably U;0,-H20, U;0,-2H2,0 and U;0,-3H,0 are formed. The 
increase in weight is related to methods of preparation, storage period 
and conditions during storage. It is shown that the process of water 
absorption is irreversible. Therefore, it is advisable to store the U,;0, 
in an ampoule, or to ignite the U;0, at 850°C in a muffle furnace for two 


hours vefore use. 
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Pharmacology 


COMPARATIVE STUDY OF PHOTODYNAMIC EFFECT, PHOTOSENSITIVE SIDE ACTION AND 
TOXICITY OF SEVERAL CHINESE PHOTOSENSITIZERS 


Beijing YAOXUE XUEBAO [ACTA PHARMACEUTICA SINICA] in Chinese Vol 21 No l, 
29 Jan 86 pp 7-11 


{English abstract of article by Wang Naiqin [3769 5082 0530], et al., of 
Beijing Institute for Cancer Research, Beijing; and Ye Ailian [0673 1947 5571], 
et al., of Beijing Pharmaceutical Industrial Institute, Beijing] 


[Text] In this paper, the photodynamic effect, photosensitive side action and 
acute toxicity of three new Chinese hematoporphyrin derivatives, HPS, LF-019 
and Y-HPD, are compared with those of America-made photofrin II. The 
experimental method included: (1) the killing effect of the photosensitizer 
plus exposure to light on the MGC-803 cells; (2) the survival time of mice 
after receiving L¢,5 cells pretreated with the photosensitizer and light; 

(3) the inhibitory effect on *H-TdR and *H-UR incorporated into DNA and RNA; 
(4) the uptake and storage of photosensitizers in L)2;9 cells; (5) the ear 

and skin indexes: showing skin photosensitive side action and (6) the LDso 

in mice. 


The results show that the photodynamic effect of the three Chinese photo- 
sensitizers waS-more potent-than that of photofrin II. In the three Chinese 
photosensitizers, except for the skin photosensitive side action, Y-HPD was 
more potent than HPS or LF-019 in the killing effect of MGC-803 cells, the 
survival time of L¢,5 mice, the inhibitory effect of nucleic acid biosynthesis 
and the content of photosensitizers in L,;2;9 cells. The acute toxicity of 
Y-HPD was lower than that of other photosensitizers. (Paper received 

4 May 1985.) 
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STUDIES OF SYNTHESIS AND STRUCTURE-ACTIVITY RELATIONSHIPS OF N-(a,8- 
DISUBSTITUTED—4-CHLOROC I NNAMYL )-SEC-BUTYLAMINES 


Beijing YAOXUE XUEBAO [ACTA PHARMACEUTICA SINICA] in Chinese Vol 21 No l, 
29 Jan 86 pp 20-28 


[English abstract of article by Peng Shiqi [1756 1597 1142], et al., of the 
College of Pharmaceutical Sciences, Beijing Medical University, Beijing; 
and Chen Suming [7115 4790 2494], et al., of the Institute of Chemistry, 
Chinese Academy of Sciences, Beijing] 


[Text] In this paper the synthesis of 13 N(a, 8-disubstituted-4-chlorocin- 
namyl)-sec-butylamine derivatives (4Ea, 4Eb~e, 4Zb~e, 6Ea~b, 6Za~b) are 
reported. The configuration of pairs of geométrical isomers has been 
assigned by NMR. Four compounds have been configurationally confirmed by 
X-ray diffraction. The 6, (the twist angle between the phenyl ring and the 
carbon-carbon double bond) and 62 (the twist angle between the carbonyl group 
and the carbon-carbon double bond) of stable conformations of the compounds 
have been calculated by the EHMO method. The values of 6,, 62, and anti- 
convulsant activity (EDs,) of 18 compounds are given here. The relationships 
of 6, and §2 and the anticonvulsant activities of these compounds are 
discussed. The pattern in which these compounds bind with the receptor in 
order to produce the anticonvulsant effect is tentatively proposed. (Paper 
received 12 February 1985.) 
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